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SECTION  I 

INTRODUCTION  AND  BACKGROUND 

This  report  addresses  the  problem  of  determining  the  conditions 
under  which  the  flow  in  a cooled  laminar  boundary  layer  becomes  unstable 
when  excited  by  low  amplitude  disturbances.  The  flow  is  considered  to  be 
unstable  when  a low  amplitude  disturbance  imbedded  in  a cooled  laminar 
boundary  layer  flow  grows  in  time  and/or  space.  The  physical  significance 
of  this  problem  is  that  such  processes  initiate  the  transition  from  laminar 
to  turbulent  flow  conditions  in  the  boundary  layer.  The  skin  friction  and 
rate  of  heat  transfer  are  typically  an  order-of-magnitude  greater  in  turbulent 
boundary  layers  than  in  laminar  boundary  layers  when  compared  at  the  same 
Reynolds  number.  Therefore,  the  ability  to  understand,  predict,  and 
possibly  suppress,  boundary  layer  transition  is  of  great  engineering  importance 
in  the  design  of  flight  vehicles. 

The  problem  of  predicting  the  conditions  under  which  a laminar 
boundary  layer  will  begin  to  undergo  the  transition  process  represents 
one  of  the  great,  and  as  yet  unsolved,  problems  of  fluid  dynamics.  The 
mathematical  description  of  fluid  dynamic  phenomena  is  a system  of  non- 
linear partial  differential  equations.  To  the  present  time,  this  system  has 
resisted  all  attempts  to  obtain  solutions  which  are  valid  for  boundary 
layer  flows  undergoing  transition.  Therefore,  simplified  versions  of  the 
general  equations  are  used  to  study  the  transition  process.  The  simplified 
equations  used  in  the  present  study  are  called  the  small  disturbance  stability 
equations.  They  are  a linearized  version  of  the  general  (Navier-Stokes) 
equations.  The  linearized  small  disturbance  equations  are  themselves 
not  very  amenable  to  direct  solution  because  of  the  partial  differential 
nature  of  the  equations.  Therefore,  a further  simplification  is  introduced. 

The  boundary  layer  is  considered  to  be  locally  parallel.  That  is,  the  growth 
of  the  boundary  layer  thickness  over  a wavelength  of  a disturbance  is 
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neglected.  This  simplification  allows  one  to  reduce  the  partial  differential 
equations  to  ordinary  differential  equations.  Thus,  the  parallel  flow  boundary 
layer  stability  equations  have  been  used  in  the  work  described  in  this  report. 
This  model  was  employed  to  study  the  conditions  under  which  small 
disturbances  will  grow  in  space  and/or  time  in  a cooled  laminar  boundary 
layer  flow . 

The  reasons  why  cooled  laminar  boundary  layers  are  specifically 
addressed  in  this  report  are: 

1.  Previous  experimental  investigations  of  boundary  layer 
transition  have  shown  that  cooling  a boundary  layer  generally  (but  not 
always)  increases  the  stability  of  the  boundary  layer.  Thus,  the 
potential  exists  for  using  boundary  layer  cooling  as  a means  for 
delaying  the  onset  of  transition. 

2.  Recent  advances  in  the  field  of  parallel  flow  boundary  layer 
stability  theory  have  made  it  possible  to  further  advance  the  use  of 
this  theory  as  a method  for  predicting  the  onset  of  transition. 

3.  Recently  published  measurements  of  transition  in  shock  induced 
flows  have  shown  an  intricate  variation  of  transition  with  wall  cooling. 
These  observations  and  other  experimental  results  cannot  be 
reasonably  explained  with  existing  theoretical  results.  New  theoretical 
results  using  the  most  powerful  methods  currently  available  for 
solving  the  stability  equations  are  necessary  to  understand  these 
observations. 

4.  The  present  energy  crisis  has  provided  a new  impetus  for 
investigating  both  new  and  old  techniques  which  could  be  used  to 
reduce  the  aerodynamic  drag  of  flight  vehicles.  Boundary  layer 
control  through  wall  cooling  thus  might  be  a useful  means  for  reducing 
drag,  conserving  fuel  resources,  and  extending  the  range  and  payload 
of  aircraft. 
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1.1  BOUNDARY  LAYER  STABILITY  TERMINOLOGY 

In  this  paragraph  a number  of  the  more  important  terms  which  are  used 
to  classify  the  different  types  of  disturbances  which  occur  in  compressible 
boundary  layer  flows  are  defined. 

1.1.1  Sonic,  Subsonic,  and  Supersonic  Disturbances 


A small  disturbance  having  a velocity  relative  to  the  free- stream 
which  is  equal  in  magnitude  to  the  speed  of  sound  in  the  free-stream  is 
called  a sonic  (acoustic)  disturbance.  The  following  definitions  are  introduced; 


c^  = velocity  of  a disturbance  (phase  velocity), 
a = speed  of  sound  in  the  free-stream. 


U = velocity  of  the  free-stream, 


M = Mach  number  of  the  free-stream. 


The  velocity  of  the  free-stream,  relative  to  an  observer 

5je 

moving  with  the  disturbance  velocity,  is  U - c . The  relative  Mach 
• ce  r 

number  (M)  is; 

M = (U  - c )/  a , (1) 

00  00  , 1*  00 

jC 

which  after  multiplying  and  dividing  by  U , becomes 

00 

A 

M = ( 1 - c ) M , (2) 

oo  r ® 

* * 

where  c is  the  nondimensional  phase  velocity  c /U  . From  the  definitions 
r r ® 

of  a sonic  disturbance  given  at  the  beginning  of  Section  1.1.1,  it  follows 

that  M = + 1 for  a sonic  disturbance.  Thus,  if  Equation  (2)  is  solved  for 

00  “ 

c^  cor  responding  to  a sonic  disturbance,  the  following  two  expressions  for 
the  phase  velocity  of  sonic  disturbances  are  obtained; 


c = 1+1/M  , M = - 1 , 
r * 

c = 1 - 1/M  , M = + 1. 

r oo 
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Thus,  for  subsonic  disturbances; 


(1  - 1/M  ) < c < ( 1 + 1/M  ) (5) 

oo  r ® 

A 

For  supersonic  disturbances  I M ^ | > 1 and  the  phase  velocity  is; 

A 

c > 1 + 1/M  , M < -1  <6> 

r ® 

A 

c < 1-1/M  , M > +1  (7) 

I*  00 

1.1.2  Disturbances  Having  a Supersonic  Relative  Flow  Regtgfl. 

Disturbances  were  classified  in  Section  1.1.1  as  subsonic, 
sonic  or  supersonic  based  on  the  speed  of  sound  in  the  free- stream. 

Because  of  the  temperature  gradient  which  exists  in  compressible  boundary 
layers,  it  is  possible  to  classify  disturbances  as  subsonic,  sonic,  or  super- 
sonic with  respect  to  the  speed  of  sound  at  any  point  in  the  boundary  layer. 
Thus,  M (n)  at  any  point  r\  in  the  boundary  layer  can  be  defined  as  the 
velocity  of  the  flow  at  r\  relative  to  the  velocity  of  the  disturbance  divided 
by  the  local  speed  of  sound,  or; 

M(n)  = (u"  (n)  - c'r  )/ a' (n).  (8) 


1.  1.3  Incoming  and  Outgoing  Waves 


A supersonic  low  amplitude  disturbance  produced  in  the 
free- stream  will  have  an  associated  set  of  Mach  waves  as  shown  in  Figure  1. 
If  a boundary  layer  is  present  in  the  flow,  Mach  waves  will  impinge  on  it. 
Such  a Mach  wave  impinging  on  a boundary  layer  is  called  an  incoming  wave. 
The  wave  produced  by  reflection  of  an  incoming  wave  on  the  boundary  is 
called  an  outgoing  wave  as  illustrated  in  Figure  1.  The  reflected  wave  can 
be  considered  to  have  the  properties  of  a Mach  wave  of  opposite  sense  to  the 
incoming  Mach  wave. 
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Figure  1.  Sketch  Showing  Defini'ion  of  Incoming  and 

Outgoing  Waves,  (a)  c <1-  l/M  , (b) 

c > 1 + 1/M  . r 

r oo 


1.1.4 


In  the  parallel  flow  linear  stability  theory,  the  fluctuation 
of  any  disturbance  parameter  E is  defined  by  a traveling  wave  of  the  form; 

■ci  „ i \ i(or,x+az-u)t)  /ii\ 

E (x  , y,  z , t)  = e (y)  e 1 3 ' (11) 

If  a 3 is  zero.  Equation  (11)  describes  a plane  wave  which  propagates  in  the 
direction  of  the  free-stream  flow.  These  waves  are  called  two-dimensional 
waves  (2-D).  If  a3  is  not  zero.  Equation  (11)  describes  a wave  which 
propagates  obliquely  to  the  free-stream  flow.  Such  a wave  is  called  a three- 
dimensional  wave  (3-D).  The  angle  between  the  direction  of  propagation 
and  the  free-stream  flow  is  called  the  wave  angle  (0  ). 


1.1.5  Critical  Reynolds  Number 

Generally  there  exists  a streamwise  location  in  a boundary 
layer  such  that  between  the  leading  edge  of  an  aerodynamic  body  and  this 
location  all  disturbances  are  damped.  The  Reynolds  number  based  on  the 
distance  to  this  location  from  the  leading  edge  is  called  the  critical  Reynolds 


number. 
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1.2  REVIEW  OF  COOLED  LAMINAR  BOUNDARY  LAY  ER  STABILITY 
. During  the  past  decade,  great  strides  have  been  made  toward  developing 

< 

a better  understanding  of  the  causes  and  nature  of  transition  in  wind  tunnel 
and  free-flight  boundary  layers.  The  greatest  analytical  progress  towards 
a better  understanding  of  transition  has  been  through  the  application  of  the 
linear  stability  theory  of  parallel  flows. 

Mack  ^ ^ made  a major  contribution  in  the  application  of  the 
linear  stability  theory  as  a tool  for  understanding  boundary  layer  transition 
by  di  scovering  multiple  mode,  multiple  families  of  solutions  to  the  compressible 
flow,  small  disturbance  equations.  He  demonstrated  the  physical  importance 
of  these  new  types  of  solutions  and  formulated  a quasi-theoretical  model 
of  transition  for  supersonic  flows  based  on  the  linear  stability  theory  of 
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parallel  flows.  Mack's  quasi- theoretical  model  of  transition  includes  two 
mechanisms  for  relating  the  influence  of  environmental  disturbances  to  the 
response  of  the  boundary  layer.  One  mechanism  is  that  of  the  forced 
oscillation  of  the  boundary  layer  caused  by  a specific  type  of  free  stream 
pressure  disturbance.  The  other  mechanism  is  the  response  of  the  boundary 
layer  to  a spectrum  of  pressure  disturbances.  Within  the  context  of  these 
two  mechanisms  he  was  able  to  study  the  effect  of  Mach  number,  unit 
Reynolds  number  (Reynolds  No.  per  foot  of  the  free- stream),  and  wall 
cooling  on  transition  in  wind  tunnel  and  free-flight  boundary  layers.  He 
was  basically  successful  in  explaining  how  environmental  disturbances  can 
affect  transition  as  Mach  number,  unit  Reynolds  number,  and  wall  cooling 
change.  The  accuracy  of  Mack's  transition  model  appears  to  be  limited  , 
for  the  most  part,  by  the  accuracy  with  which  the  nature  of  the  environmental 
disturbances  themselves  are  known.  In  order  to  apply  Mack's  model  to  the 
problem  of  predicting  the  location  of  the  onset  of  transition  in  a flow 
situation,  one  must  have  available  the  basic  results  of  a parallel  flow, 
small  disturbance  stability  analysis;  namely,  neutral  stability  curves  and 
amplification  rate  versus  frequency  for  a range  of  Reynolds  numbers. 

The  considerable  progress  that  has  been  made  toward  understanding 
transition  in  wind-tunnel  and  free-flight  boundary  layers  has  not  been 
matched  by  progress  in  understanding  transition  in  laminar  boundary 
layers  induced  by  a moving  shock  wave.  This  situation  persists  in  spite  of 
the  fact  that  a considerable  amount  of  experimental  data  on  transition  in 

shock-tube  boundary  layers  has  been  generated  during  the  past  two  decades. 

(2) 

Morkovin  has  summarized  the  nuances  and  contradictory  observations 

that  have  been  made  concerning  transition  in  shock  tube  facilities  and 

has  also  conducted  an  assessment  of  the  level  of  understanding  of  wind- 

. . (3) 

tunnel  and  free-flight  boundary  layer  transition 

(4) 

Recently,  Boison  conducted  a series  of  careful  experiments  which 
considered  shock  tube  boundary  layer  transition.  He  examined  the  influence 
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on  transition  of  environmental  factors,  such  as  wall  vibration  and  free 

stream  acoustic  disturbances,  in  addition  to  the  influence  of  wall  cooling  , 

unit  Reynolds  number,  and  Mach  number.  In  these  experiments,  extraordinary 

precautions  were  taken  to  eliminate  wall  roughness  as  a possible  cause  of 

early  transition.  The  transition  is  usually  expressed  in  terms  of  the 

"transition"  Reynolds  number  (R^  ) ^ based  on  the  length  of  the  boundary 

layer  at  transition.  Boison's  data  at  a free  stream  Reynolds  number  per 

foot  of  5 x 105  showed  that,  as  wall  cooling  is  increased,  the  transition 

Reynolds  number  first  increased,  then  decreased  (reversal),  then  increased 

again  (re- reversal) , then  decreased,  increased  again,  and  so  on;  showing  a 

number  of  transition  reversals  and  re- reversals  . The  transition  Reynolds 

number  was  found  to  be  extremely  sensitive  to  small  changes  in  cooling 

especially  at  the  larger  cooling  rates.  Figure  2 shows  Boison’s  transition 

5 

results  for  R /ft  = 5 x 10  along  with  results  of  several  other  investigators. 

e 

The  wall  cooling  rate  is  expressed  in  terms  of  the  wall  temperature  T^ 

and  the  boundary  layer  edge  temperature  Tg  in  this  figure.  The  ratio 

T / T is  inversely  related  to  the  magnitude  of  wall  cooling  (decreasing 

T / T corresponds  to  increasing  wall  cooling), 
w e 

The  first  transition  reversal  loop  (0.17  < T /T  0.36)  seems  to  be 

r we 

fairly  well  established  although  more  results  would  be  desirable  to 
definitely  show  that  only  one  loop,  rather  than  perhaps  two  loops,  exists 
in  this  wall  cooling  range.  Also  shown  in  Figure  2 is  the  calculation 

(5) 

of  Reshotko  for  the  wall  cooling  required  for  complete  stabilization  of  2-D 
disturbances.  It  should  be  noted  that  this  result  is  not  based  on  a 
completely  valid  solution  of  the  stability  equations  as  determined  by  Mack 
(see  page  14.1  of  Reference  6).  In  particular,  mean  velocity  and 
temperature  profiles  characteristic  of  shock  tube  boundary  layers  were 
not  used. 
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Figure  2.  Transition  Reynolds  No.  Correlation  for  Re /ft  = 5 x 10  . 

The  description  of  only  one  stability  analysis,  using  velocity  and 

temperature  profiles  characteristic  of  compressible  boundary  layers 

behind  a moving  shock  appears  in  the  literature.  That  analysis,  performed 

(7) 

by  Ostrach  and  Thornton  in  1962,  produced  minimum  critical  Reynolds 
numbers  that  greatly  exceeded  the  measured  transition  Reynolds  numbers  - 
the  deviation  increasing  by  orders  of  magnitude  with  increasing  cooling,  i.  e., 
higher  shock  Mach  numbers.  Their  analysis  was  based  on  the  Dunn-Lin 
asymptotic  approach  ' ' and  did  not  take  into  account  3-D  disturbances. 
Ostrach  and  Thornton  concluded  that  shock  tube  transition  did  not  occur 
through  amplification  of  low  amplitude  waves  but  rather  was  caused  by  large 
environmental  disturbances  which  are  present  in  shock  tubes. 

The  existence  of  multiple  mode,  multiple  family  solutions  to  the 

(2) 

parallel  flow  linear  stability  equations  led  to  the  suggestion  by  Morkovin 
that  the  Ostrach -Thornton  analysis  be  redone  using  the  full  set  of  parallel 
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flow,  linear  stability  equations.  Mack  had  shown  that,  depending  on 

the  amount  of  wall  cooling  (T  / T ),  2-D  first  mode  unstable  disturbances 

w e 

could  be  stabilized  and  second  mode  2-D  disturbances  destabilized. 

Morkovin  ^ thought  that  perhaps  the  intricate  cooling  effects  on  transition 

in  shock  tube  boundary  layers  could  be  explained  as  the  stabiliziation  of 

first  mode  disturbances  over  certain  ranges  of  T /T  and  destabilization  of 

we 

higher  modes  over  other  ranges  of  T /T  . This  suggestion  was  put  forth 

W 6 

prior  to  Boison's  work.  When  Boison  discovered  the  existence  of  the 
multiple  transition  reversal  loops  in  his  shock  tube  transition  data  he  saw 
the  possibility  that  these  loops  might  be  explained  through  a stability 
analysis  of  the  shock  tube  laminar  boundary  layer  which  recognized  the 
existence  of  multiple  mode  solutions. 

In  early  1972,  an  effort  was  begun  by  one  of  the  authors  (Boehman)  to 
reconduct  the  Ostrach- Thornton  analysis  using  the  authors'  computer 
programs  which  had  been  successfully  applied  to  the  steady  flow  supersonic 
boundary  layer  stability  problem  12^.  This  initial  effort  was  disappointing 
insofar  as  higher  mode  solutions  for  the  shock  tube  laminar  boundary  layer 
could  not  be  found.  The  first  mode  2-D  solutions  only  confirmed  the  results 
of  Ostrach  and  Thornton.  In  this  initial  effort,  it  was  found  that  the  higher 
modes  are  associated  with  large  values  of  the  product  of  wave  number 
(reciprocal  of  the  spatial  wavelength)  and  Reynolds  number.  This  is 
characteristic  of  higher  modes  in  low  free-stream  Mach  number  flows  with 
high  cooling  rates.  It  was  also  found  that  larger  numbers  of  very  small 
integration  steps  were  required  to  integrate  the  stability  equations  for  shock- 
tube  boundary  layers.  This  report  describes  subsequent  efforts  to  locate 
higher  mode  solutions  to  the  parallel  flow  stability  equations  for  highly  cooled 
compressible  boundary  layers.  j: 


1.3  REPORT  OUTLINE 

Section  II  of  this  report  describes  the  analysis  of  stability  in  shock 
induced  boundary  layer  flows.  The  section  begins  with  a description  of  the 
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calculation  of  the  boundary  layer  mean  flow  profiles.  This  description 
is  followed  by  a treatment  of  the  formulation  of  the  stability  equations. 
Finally,  some  of  the  limits  of  applicability  of  the  formulation  are  reviewed. 

Section  III  presents  the  results  of  the  analysis  of  boundary  layer 
stability  in  shock  induced  flows.  The  results  are  preceeded  by  a general 
description  of  the  families  of  solutions  identified  by  Mack.  The  specific 
results  of  this  investigation  are  then  presented  and  discussed. 

f 1 

Section  IV  describes  the  application  of  parallel  boundary  layer 
stability  theory  to  another  class  of  cooled  boundary  layer  flows  - steady 
flow  boundary  layers.  A short  description  of  the  problem  formulation 
is  followed  by  a presentation  and  discussion  of  the  results.  The  important 
similarities  and  differences  with  the  results  from  shock  induced  flows  are 
then  discussed. 

In  Section  V,  the  important  conclusions  of  the  study  are  presented  for 
j both  shock  induced  and  steady  flow  boundary  layers  and  recommendations 

for  further  work  are  presented. 
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SECTION  II 


PARALLEL  STEADY  FLOW  ANALYSIS 
OF  SHOCK  INDUCED  FLOWS 

The  methods  used  to  investigate  the  stability  of  the  shock  tube  laminar 
boundary  layers  are  presented  in  this  section.  In  Section  2.1  the  methods 
used  to  generate  the  mean  flow  velocity  and  temperature  profiles  are  presented. 
In  Section  2.2  it  is  shown  that,  when  quasi-parallel  , quasi -steady  approxi- 
mations are  made  for  shock  induced  boundary  layers,  a set  of  stability 
equations  are  obtained  which  are  identical  to  the  stability  equations  for 
steady  boundary  layer  flows. 

A computer  program,  previously  developed  for  steady  flow  boundary 
layer  stability  analysis,  can  then  be  applied  to  shock  induced  boundary 
layers.  The  equations  and  solution  techniques  of  this  program  are  described 
in  detail  elsewhere  , but  a brief  description  is  included  in  Appendix 

B.  Section  2.3  discusses  the  range  of  flow  parameters  that  were  considered 
in  this  investigation. 

2.1  MEAN  FLOW  PROFILES 


Mean  flow  velocity  and  temperature  profiles  are  obtained  by  numerically 
integrating  the  laminar,  compressible  flow,  boundary  layer  equations 
written  in  a shock  fixed  coordinate  system.  In  this  coordinate  system 
described  in  Figure  3 the  flow  is  steady  and  the  boundary  layer  equations 

(13) 

are  formulated  and  solved  according  to  the  procedures  developed  by  Mack 
which  served  as  the  basis  of  the  authors'  computer  program.  Specifically, 
the  program  is  set  up  to  solve  the  steady  laminar,  compressible  flow, 
boundary  layer  equations  for  flat  plates  with  and  without  heat  transfer  at 
the  wall.  For  the  shock  tube  boundary  layer  calculation,  the  input 
quantities  are  as  follows: 


• Shock  Mach  number  (M  ), 

s . 

• Undisturbed  flow  temperature  (T 


* 

) and  pressure  (p ^ 
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(o)  STEADY  COORDINATE  SYSTEM  (SHOCK  FIXED) 


(b)  UNSTEADY  COORDINATE  SYSTEM  (LAB  FIXED) 


Figure  3.  Boundary  Layer  Mean  Flow  Velocity  Profile. 

• Initial  guesses  for  the  slope  of  the  velocity  profile  at  the 
wall  and  the  heat  transfer  rate  at  the  wall. 

The  normal  shock  equations  are  solved  to  determine  the  free  stream 

>\c  ;Js  5|c  ;|c 

conditions  of  the  flow  downstream  from  the  shockwave:  M , U ,T  ,p  ,h 

e e e e e 

(enthalpy)  and  h^  (stagnation  enthalpy  of  the  free  stream  in  the  shock  fixed 
coordinate  system).  These  quantities  are  used  to  nondimensionalize 
the  boundary  layer  equations.  A basic  assumption  made  in  solving  the 
boundary  layer  equations  is  that  the  initial  temperature  of  the  wall  (T  ) 

■f  sis  W 

is  the  same  as  T ^ and  that  T^  remains  constant.  This  assumption  is 
reasonable  because  the  time  for  heat  transfer  is  very  small  (milliseconds 
or  less)  even  though  the  heat  transfer  rate  is  high  . The  program  is 
presently  set  up  for  air  which  is  considered  to  be  a gas  that  obeys  the 
perfect  gas  law  and  has  variable  specific  heats  and  transport  properties. 


13 


i I 


Because  dissociation  and  real  gas  effects  are  not  considered,  the  program 

is  only  valid  for  M < 6.  5. 

s 

Once  the  free  stream  conditions  of  the  flow  downstream  from  the 

normal  shock  have  been  determined,  the  boundary  layer  equations  are 

integrated  and  the  mean  flow  profiles  are  determined.  Some  representative 

nondimensionalized  velocity  profiles  (in  the  shock-fixed  coordinate  system) 

are  presented  in  Figure  4.  Temperature  profiles  are  presented  in  Figure  5. 

The  profiles  are  presented  in  terms  of  the  Blasius  similarity  variable 
/ Y /Ug  X . l/2 \ 

rn  - ^ ^ — J Iwhich  is  the  "nondimensionalized"  boundary  layer 

thickness.  e 

The  basic  unknowns  to  be  determined  when  solving  the  boundary  layer 
equations,  in  addition  to  the  profiles  themselves,  are  the  slope  of  the 
velocity  profile  at  the  wall  and  the  heat  transfer  rate  at  the  wall.  The 
Newton-Raphson  technique  is  used  to  search  for  values  of  these  two 
unknowns  which  will  satisfy  the  proper  boundary  conditions  at  the  edge  of 
the  boundary  layer.  After  the  solution  converges,  the  velocity  profile 
is  transformed  to  the  laboratory-fixed  coordinate  system  and  the  velocity 
profiles  in  this  coordinate  system  and  the  temperature  profiles  are  then 
punched  onto  cards  for  use  as  a data  deck  for  the  Boundary  Layer  Stability 
computer  program  (BLSTAB). 


...  ' 


- \ 


2.  2 FORMULATION  OF  THE  STABILITY  EQUATIONS  FOR 
A SHOCK  INDUCED  BOUNDARY  LAYER 

The  flow  in  a shock  induced  boundary  layer  is  extremely  difficult  to 
analyze.  Simplifying  assumptions  concerning  the  flow  must  be  made  in 
order  to  formulate  the  problem  in  a practically  solvable  form.  There  are  a 
variety  of  methods  of  formulating  the  problem,  but  care  must  be  taken  to 
insure  that  a consistent  set  of  simplifying  assumptions  and  nondimensionali- 
zations  are  made. 
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The  stability  equations  for  shock  induced  boundary  layer  flow  were 

(7) 

formulated  by  Ostrach  and  Thornton  . They  employed  mean  flow 

(14) 

boundary  layer  profiles  as  calculated  by  Mirels  in  a shock  fixed 
coordinate  system.  Ostrach  and  Thornton  then  wrote  the  Dunn- Lin  stability 
equations  for  a coordinate  system  fixed  on  the  moving  shock  wave  and 
normalized  all  velocities  (mean  flow  profile,  disturbance  velocity  components, 
and  phase  velocity)  by  the  velocity  . This  procedure  produced  a set 
of  stability  equations  that  were  identical  to  the  Dunn- Lin  equations**  for 
steady  flow  except  that  Ostrach  and  Thornton  "redefined"  Mach  and  Reynolds 
numbers.  The  redefined  free  stream  Mach  number  was  considered  to  be  the 
Mach  number  of  the  shocked  gas  in  the  laboratory  fixed  coordinate  system 
(M^)  rather  than  the  corresponding  Mach  number  in  the  shock-fixed  coordinate 
system  (M^).  The  redefined  Reynolds  number  was  based  upon:  (1)  the  free 
stream  velocity  of  the  gas  in  the  laboratory-fixed  coordinate  system  (U^) 
instead  of  the  velocity  in  the  shock  fixed  coordinate  system  (Ue  ),  and  (2) 
a characteristic  length  equal  to  the  distance  traveled  by  a free  stream  fluid 

element  in  the  laboratory  fixed  coordinate  system  (x  ) during  the  total  time 

* 

that  the  fluid  element  is  in  motion  (t  ),  instead  of  distance  that  the  fluid 

* 

element  is  behind  the  shock  wave  (X  ).  Since  these  choices  are  not  consistent, 
the  subsequent  formulation  was  examined  in  some  detail. 

The  formulation  of  the  stability  equations  by  Ostrach  and  Thornton 
is  not  straightforward.  A critique  of  their  formulation  is  presented  in 
Appendix  A.  The  essence  of  their  approach,  however,  is  contained  in  the 
following  statements. 

1.  The  nondimensional  Dunn-Lin  version  of  the  small  disturbance 
equations  were  written  for  a shock  fixed  coordinate  system. 

2.  The  velocities  in  these  equations  were  normalized  with  respect 
to  the  free  stream  velocity  in  a laboratory  fixed  coordinate  system 


The  Dunn-Lin  equations  are  not  a complete  set  of  the  small  disturbance 
stability  equations.  Certain  terms  are  omitted  to  permit  3-D  solutions  to  be 
obtained  from  2-D  solutions. 


I 
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# $ * * 

(U  , where  U_  = U - U ) rather  than  by  the  free  stream  velocity  in 
c.  4 s e . 

5jC 

the  shock  fixed  coordinate  system  (U  ). 

sjc 

3.  The  characteristic  length  ( l \ ) was  taken  to  be  proportional  to 

the  local  boundary  layer  thickness  in  the  shock  fixed  coordinate  system, 

* * .M2 

given  by:  a - ( 2 X v*  / U',s)  1 . 

w e 

yV~  ^ ^ ^ 

4.  The  characteristic  time  (t  ) was  given  by,  t = l /U  - U ). 

c c s e 

5.  Fluid,  thermodynamic,  and  transport  properties  were  non- 
dimensionalized  with  their  corresponding  free  stream  values. 

When  the  dimensional  Dunn- Lin  small  disturbance  equations  are  non- 
dimensionalized  using  the  characteristic  quantities  defined  above,  the  final 
form  of  the  Ostrach  and  Thornton  stability  equations  is  obtained  (Equations 
9a  through  13a  and  Equation  14  in  Reference  7). 

The  nondimensionalizations  and  normalizations  described  above 

effectively  transformed  the  stability  equation  formulation  out  of  the  shock 

fixed  coordinate  system  into  the  laboratory  fixed  coordinate  system. 

(14) 

However,  the  mean  flow  profiles  as  derived  from  Mirels1  analysis 
remained  in  the  shock  fixed  coordinate  system.  The  stability  formulation 
was  now  inconsistent  with  the  mean  flow  formulation.  This  inconsistency 
was  resolved  by  "redefining"  certain  nondimensionalizations  (such  as 
Reynolds  number).  This  "redefinition"  effectively  transformed  the  stability 
analysis  back  into  shock  fixed  coordinates.  The  final  formulation  was 
therefore,  entirely  in  shock  fixed  coordinates.  A more  detailed  critique 
of  this  formulation  is  contained  in  Appendix  A. 

The  final  formulation  of  the  stability  equations  by  Ostrach  and 
Thornton  was  unsuitable  for  use  with  existing  computational  programs  at 
the  University  of  Dayton  which  were  constructed  in  a laboratory  fixed 
coordinate  system.  Therefore,  a different  approach  to  formulating  the 
stability  equations  was  developed.  This  approach  starts  with  the 


rr 


■ 
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linearized  complete  small  disturbance  equations  in  dimensional  form  (rather 
than  the  Dunn- Lin  equations)  written  in  a laboratory  fixed  coordinate  system. 
Then  both  the  parallel  flow  and  quasi- steady  approximations  were  invoked. 
Finally,  all  dependent  variables  were  nondimensionalized  with  respect  to 
free  stream  quantities  and  all  spatial  variables  were  nondimensionalized 
with  respect  to  a scale  length  proportional  to  the  boundary  layer  thickness. 
This  procedure  yields  a system  of  stability  equations  which  are  identical 


to  the  Ostrach  and  Thornton  stability  equations  except  that  all  parallel  flow 

1/2 

i 

r " ' ■ x 


terms  are  included  and  the  Reynolds  number  (R  ) is  given  by  (R  ) 

1/2  * * 1 X 
instead  of  (2  R ) (T  / T ).  (The  Ostrach-Thornton  formulation  includes 


1/2 

x we 

only  those  terms  contained  in  the  Dunn-Lin  stability  equations. 


In  summary,  the  stability  equations  used  in  the  present  analysis  are 
the  same  as  the  stability  equations  for  steady  flow  over  a flat  plate. 

The  stability  analysis  is  thus  a stability  analysis  of  a given  mean  flow  profile 
and  is  valid  over  changes  in  distances  and  times  that  are  small  enough 
so  that  the  change  of  profiles  in  both  time  and  space  can  be  ignored. 

The  computer  program  which  implements  this  formulation  is  described  in 
Appendix  B. 


2.  3 


RANGE  OF  FLOWS  ANALYZED 


Eigensolutions  to  the  linear  stability  equations,  with  the  parallel 
flow  assumption,  were  computed  for  both  a quasi- steady  flow  shock  tube 
system,  and  a steady  flow  wind  tunnel  system.  The  effects  of  wall  cooling 
on  boundary  layer  stability  in  shock  tube  type  flows  were  examined  for  both 


subsonic  (Mg  =1.5,  M^  = 0.  6027  and  supersonic  (Mg  s 5.  0,  < I.  841 ) flows. 


The  effects  of  wall  cooling  on  boundary  layer  stability  for  Mach  numbers 
less  than  0.9  in  steady  wind  tunnel  type  flows  were  also  investigated. 


2.  3.  1 


Shock  Tube  Flows 


The  shock  Mach  numbers  considered  in  this  analysis  extended 


from  M = 2.  50  to  M = 5.00  in  increments  of  0.25.  In  addition,  treatment 
s s 


i| 

i 

<1 
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was  given  to  M = 1.5,  as  previously  noted.  These  calculations  were 
s 

performed  for  the  "standard"  wall  cooling  ratios  (T  /Te),  as  given  in 

Table  1.  The  "standard"  wall  cooling  is  that  which  occurs  when  the  shock 

tube  walls  remain  at  ambient  temperature  regardless  of  the  temperature 

of  the  shocked  gas.  In  addition,  calculations  were  performed  at  shock  Mach 

numbers  of  M = 2.  5 and  M =2.75  with  specified  "nonstandard"  amounts 
s s 

of  wall  cooling.  Numerical  solutions  to  all  the  above  situations  were 

generated  for  the  case  of  two-dimensional  disturbances.  Three  dimensional, 

oblique  solutions  were  obtained  for  only  the  cases  of  M^  =1.5  (subsonic) 

and  M =3.0  (supersonic)  flows.  The  results  of  these  calculations  are 
s 

presented  in  Section  III. 

It  should  be  noted  that  the  technique  employed  is  valid  only  for 
shock  tube  flows  where,  in  general,  M < 6.  5,  as  "real  gas"  effects  have  been 
ignored. 

2.  3.  2 Subsonic  Wind  Tunnel  Flows 

Two  and  three  dimensional  disturbances  in  steady  flow,  wind 

tunnel  type  boundary  layers  were  investigated  for  free  stream  Mach  numbers 

of  0.  603,  0.  80,  and  0.  90.  Surface  cooling  at  each  Mach  number  was  varied 

from  the  insulated  value  of  T /T  (specified  by  the  flow  conditions),  to 

w e 

T / T = 0.  620.  The  streamwise  wave  angle  of  three  dimensional  disturbances 
w e 

was  varied  from  0 to  87.  5 for  each  case.  Results  of  these  calculations 
are  presented  in  Section  IV. 
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SECTION  III 


RESULTS  OF  ANALYSIS  OF  SHOCK  INDUCED  FLOWS 

As  noted  in  Section  1.2,  Boehman  made  his  first  attempt  to  deter- 
mine higher  mode  solutions  to  the  parallel  flow  stability  equations  for  shock 
induced  flows  in  1972.  These  initial  attempts  were  unsuccessful  but  higher 
mode  solutions  were  obtained  during  the  course  of  the  present  investigation. 
The  results  of  the  search  for  higher  mode  solutions  are  presented  in  this 
section.  However,  these  higher  mode  solutionsylo  not  correlate  with 
observed  transition  phenomena.  It  was  also  pointed  out  in  Section  1.  2 
that  the  Ostrach-Thornton  investigation  of  the  stability  of  shock  induced 
flow  did  not  include  an  investigation  of  3-D  (oblique)  disturbances.  The 
search  for  3-D  unstable  disturbances  was  included  in  the  present  effort 
but  it  was  found  that  three  dimensional  disturbances  were  always  more 
stable  than  two  dimensional  disturbances  over  the  range  of  flow  parameters 
considered. 

The  principal  objective  of  the  present  effort  was  to  determine  higher 
mode  and  three  dimensional  solutions  to  the  stability  equations  for  shock 
induced  flows  in  the  expectation  that  these  solutions  would  correlate  with 
experimental  observations  of  transition  in  shock  induced  flow.  When  it 
was  found  that  this  expectation  was  not  to  be  realized,  a search  for  other 
types  of  solutions  to  the  stability  equations  was  undertaken.  A new  type  of 
solution  was  found  for  a supersonic  incoming  wave  disturbance  and  quali- 
tative correlation  with  Boison's  measurements  of  transition;  Reynolds 
number  has  been  demonstrated.  Before  proceeding  with  the  presentation 
and  discussion  of  these  new  results,  the  flow  conditions  under  which 
different  types  of  disturbances  can  be  expected  are  described. 


3.1  SUBSONIC,  SONIC,  AND  SUPERSONIC  DISTURBANCES  IN 
SHOCK  INDUCED  FLOWS 

Three  curves  of  phase  velocity  versus  shock  Mach  number  (M  ) 

A S 

are  shown  in  Figure  6.  The  curve  labeled. c (M  =-l), shows  the  minimum 

r w 
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NON-DIMENSIONAL  PHASE  VELOCITY~Cr 


M = (u-Cr)l 


Mw=-CrM2 


*y  tw 


.Cr=|+^  ( M 2 = ~ I ) 


. SUPERSONIC 


SUBSONIC  DISTURBANCES  + SuB 
HAVING  SUPERSONIC  RELATIVE 
FLOW  REGION  ( |AJ>I ) IN  BOUNDARY 
LAYER 


SUBSONIC 


fo>=o_ 
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RELATIVE  FLOW 
REGION  POSSIBLE 
4 ~ BELOW  THIS  CURVE 
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Figure  6.  Phase  Velocity  Versus  Shock  Mach  No.  in  Shock  Tube  Flows, 


phase  velocity  that  a disturbance  must  have  in  order  to  move  supersonically 

relative  to  the  flow  at  the  wall.  The  second  curve,  labeled  c = 1 + 1/M„, 

r 2 

shows  the  phase  velocity  that  a disturbance  must  have  to  move  sonically 
with  respect  to  the  disturbance  (i.  e.  , the  free  stream  is  moving  in  the 
upstream  direction  at  M - 1 relative  to  an  observer  moving  at  the  velocity 
1 + l/M^).  The  third  curve,  labeled  c = 1 - l/M^,  shows  the  phase 
velocity  that  a disturbance  must  have  to  allow  the  free  stream  to  move 
sonically  with  respect  to  disturbance  but  with  the  relative  velocity  between 
the  disturbance  and  the  free  stream  being  in  the  direction  of  flow.  These 
three  curves  together  form  a region  in  which  subsonic  disturbances  exist 
and  the  local  Mach  number  of  the  flow  at  the  wall  is  supersonic  relative  to  the 
phase  velocity.  This  region  of  phase  velocities  is  the  region  in  which 
subsonic  higher  mode  disturbances  can  exist.  Above  M^  = 3.  5,  both  subsonic 
and  supersonic  multiple  mode  solutions  can  exist. 

3.  2 INCOMING  AND  OUTGOING  DISTURBANCES 

It  was  noted  in  Section  1.  1.  3 and  illustrated  in  Figure  1 that  a pair 
of  Mach  waves  are  produced  whenever  a low  amplitude  supersonic  disturbance 
exists  in  the  free  stream.  Mack  has  noted  that  the  two  solutions  to  the 


inviscid  stability  equations  for  the  flow  in  the  free  stream  represent  a pair 
of  Mach  waves  Thus,  the  stability  equations  produce  solutions  for 

supersonic  disturbances  which  have  a well-defined  physical  interpretation. 

One  of  these  solutions  represents  a reflected  or  outgoing  Mach  wave  and  the 
other  represents  an  incident  or  incoming  Mach  wave.  In  Reference  6 Mack 
shows  that  for  neutral  outgoing  Mach  waves,  energy  is  transported  in  the 
direction  of  increasing  y (i.  e.  , into  the  free  stream)  whereas  for  neutral 
incoming  waves,  energy  is  transported  in  the  direction  of  decreasing  y, 

(i.  e. , into  the  boundary  layer). 

' 9 

For  subsonic  disturbances,  the  inviscid  free  stream  solutions  no 
longer  represent  Mach  waves,  but  instead  represent  two  exponentially 
varying  pressure  fields  caused  by  inviscid  flow  over  a moving  wavy  wall. 

Thus  for  subsonic  phase  velocities,  as  well  as  for  supersonic  phase  velocities. 
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the  solutions  to  the  inviscid  stability  equations  have  simple  physical 
interpretations.  A basic  understanding  of  the  mathematical  nature  of  the 
solutions  to  the  inviscid  stability  equations  is  helpful  for  understanding 
the  discussion  of  results  presented  in  this  section.  The  solution  to  the 

inviscid  stability  equations  for  the  free  stream  region  of  flow  has  the  general 
form  (6,: 

tt  (y ) = AeX  1 T1+  Be  X2  ^ , (12) 

where  tt  represents  the  amplitude  of  the  pressure  oscillation,  \ and  \ 
are  the  two  characteristic  values  of  the  second  order  ordinary  differential 
equation  which  describes  the  behavior  of  n in  the  free  stream  when  no 
viscous  terms  are  included  in  the  stability  equations,  r)  is  the  distance 
from  the  surface  (Blasius  similarity  variable)  and  A and  B are  constants 
of  integration.  \ ^ and  \ ^ are  represented  by 

x i = " a(l-  M22)  1/2  , (13) 


X 2 “ X 1 * (14! 

The  fundamental  difference  in  nature  between  subsonic  and  supersonic 
disturbances  is  described  in  the  following  sections  in  which  the  properties 
of  the  characteristic  values  of  neutral  disturbances  (c.  =0)  for  these  two 
cases  are  examined. 


3.  2.  1 Supersonic  Disturbances  ( | M.,  | > 1 ) 

a 2 

For  |M2|>  1,  1 - M2  is  a negative  number  and  (1-M2)^2 
a pure  imaginary  number.  Thus  \ ^ is  given  by: 

A 2 1/2 

Xj  = ± ia(M2  - 1)  (15) 


In  the  coordinate  system  moving  with  the  phase  velocity, 
the  pressure  fluctuation  p'  is  given  by 

2 a 

p*  = ioyM2(l  - cr  ) exp [ io{x  T (M2  - IP  2(T)-  T^)] ) 
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Thus,  in  the  x-p  plane,  the  above  equation  will  have  a line  of  constant 

* 

amplitude  (constant  phase)  which  is  a Mach  line  of  the  flow  relative  to  the 
phase  velocity.  The  amplitude  of  the  pressure  fluctuation  is  independent  of  r|. 
If  c > 1 + l/M^,  then  the  incoming  wave  will  correspond  to: 

~ 2 l/2 

\ = - ia  (M_  - 1 ) (incoming  wave,  c >1+  1/M_). 

12  r 2 

If  c <1  - l/M^,  the  incoming  wave  will  correspond  to: 

A 2 l/2 

\ , = + ia  (M  - 1)  (incoming  wave,  c < 1 - 1/M0). 
12  r 2 


3.2.2  Subsonic  Disturbances  (-1  < < 1) 


* 2 


1/2 

In  the  case  of  subsonic  disturbances,  (1  - M^)  is  a real 
number;  \ is  negative  and  \ is  positive.  The  pressure  fluctuation  p1  is 

1 u 

given  by: 


p'  = ioyM^(  1 - c ) exp  (icrx)  exp  |>a(l-  M2)^2(71-  T|  )]  (17) 

The  solution  whose  characteristic  value  has  a negative  real 
part  represents  the  pressure  field  over  the  moving  wavy  wall  (outgoing 
solution).  The  solution  whose  characteristic  value  has  a positive  real 
part  represents  the  pressure  field  under  the  moving  wavy  wall  (incoming 
solution). 


In  the  viscous  theory,  two  of  the  independent  solutions  are 
almost  identical  to  the  inviscid  solutions  except  for  a small  viscous  decay 
term.  Ordinarily  in  stability  theory,  the  inviscid  solution  (or  the  viscous 
counterpart)  whose  characteristic  value  has  a positive  real  part  (incoming 
solution)  is  not  used  since  this  solution  increases  exponentially  upwards 
from  the  boundary  layer  and  thus  does  not  satisfy  the  boundary  conditions  at 
infinity.  In  the  following  paragraphs  it  will  be  convenient  to  classify  the 
results  of  the  stability  investigations  in  terms  of  incoming  or  outgoing 
subsonic  or  supersonic  disturbances.  Sonic  disturbances  have  a singular 
character  as  explained  in  the  next  section. 


j 

i 
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3.  3 FAMILIES  OF  SOLUTION 


Perhaps  one  of  the  most  important  discoveries  made  by  Mack  was 

that  two  families  of  solutions  can  be  obtained  from  the  stability  equations 

for  compressible  laminar  boundary  layer  flows Mack  characterized 

the  two  families  on  the  basis  of  the  behavior  of  the  solutions  at  a =0. 

He  called  a set  of  solutions  members  of  the  c family  of  solutions  if  the 

o 

phase  velocity  approached  a value  of  c = 1 - 1/M^as  a.  approached  zero. 

He  called  a set  of  solutions  members  of  the  "regular"  family  of  solution  when 
the  phase  velocity  approached  a value  of  c^  = 1 41/M^as  0.  approached  zero. 

Most  of  the  stability  results  which  are  available  in  the  literature  for 
supersonic  flows  are  subsonic  solutions  belonging  to  the  c^  family  of 
solutions.  If  the  free  stream  Mach  number  is  less  than  one  (subsonic  boundary 
layer  flows),  the  c family  of  solutions  has  its  origin  at  c^  = 0.  Relatively 
few  "regular"  family  solutions  are  available  in  the  literature  --those  which 
do  exist  were  determined  by  Mack  or  Boehman. 

From  the  preceeding  discussion  it  is  evident  that  the  two  sonic  phase 
velocities  c = 1 - 1/M  and  c =1  +1/M  define  two  separate  families  of 
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3.4  RESULTS  FOR  SHOCK  INDUCED  FLOWS 

This  section  begins  with  a presentation  of  the  results  which  were 
obtained  from  the  first  effort  to  reconduct  the  Ostrach-Thornton  analysis. 

As  was  mentioned  in  Section  1.  1,  the  effort  was  carried  out  before  the 
present  effort  was  undertaken  but  is  included  here  for  the  sake  of  completeness. 

3.4.1  Subsonic  Outgoing  Disturbances 

As  was  mentioned  in  the  Introduction,  the  first  effort  showed 
that  very  small  integration  step  sizes  are  required  to  numerically  integrate 
the  stability  equations  for  the  shock-tube  boundary  layer  when  compared 
to  step- size  requirements  for  treating  wind  tunnel  boundary  layers.  Figure  7 
shows  the  neutral  stability  curve  for  2-D  disturbances  which  was  obtained 
for  the  M -1.5  boundary  layer.  On.-  hundr.  d and  thirty-six  integration 


1 


26 


DIMENSIONLESS  WAVE  NUMBER 


Figure  7.  Two  Dimensional  Disturbance  Neutral  Stability  Curve 
at  Mg  = 1.5  for  Shock  Tube  Boundary  Layer. 


steps  were  required  to  obtain  these  results.  The  minimum  critical  Reynolds 

(R  . ) is  seen  to  be  13,  180  and  compares  to  the  Ostrach-Thornton  value  of 
crit 

14,000  for  the  same  value  of  M . These  calculations  as  well  as  those  of 

s 

those  of  Ostrach  and  Thornton  were  for  2-D  waves.  A value  of  R .of 

8 Crit 

13,  180  corresponds  to  a value  of  R^  of  1.  737  x 10  which  is  about  two  orders 
of  magnitude  higher  than  measured  transition  Reynolds  numbers  for  = 1.5. 

The  remainder  of  the  results  presented  in  this  report  were  obtained  as  part 
of  the  current  effort. 

Since  R . for  2-D  disturbances  at  M =1.5  obviously 
crit  s 

did  not  correlate  with  measured  transition  Reynolds  numbers,  3-D 

disturbances  were  investigated  at  = 1.  5.  The  results  are  shown  in 

Figures  8 and  9 for  several  Reynolds  numbers.  These  results  show  that 

at  M =1.  5,  oblique  disturbances  are  always  more  highly  damped  than 

s - 

2-D  disturbances.  Similar  calculations  performed  for  M = 3.  0 also  show 

s 

27 


afj 


i 

i 


I 

..  ! 
I 


I 

. 


that  3-D  disturbances  are  more  highly  damped  than  2-D  disturbances. 

Typical  results  at  M =3.0  are  shown  in  Figures  10  and  11.  The  results 

shown  in  Figures  8-11  are  for  disturbances  belonging  to  the  cq  family. 

No  higher  mode  amplified  solutions  for  this  family  of  solutions  are  possible 

for  shock  tube  boundary  layers  since  the  wall  cooling  is  sufficient  to  remove 

(6) 

the  generalized  inflection  point  from  these  boundary  layers 

First  and  higher-mode  outgoing  subsonic  solutions  belonging 
to  the  "regular"  family  of  solutions,  which  are  always  amplified  when  c^  1, 
have  been  investigated.  The  regular  family  of  solutions  persist  to  low  sub- 
sonic Mach  numbers,  since  it  is  always  possible  to  find  a c^  large  enough 
so  that  a supersonic  relative  flow  region  exists  somewhere  in  the  boundary 
layer,  except  at  M = 0.  In  Figure  6,  the  curve  labeled  c^  = 1 + 1 /M^ 

A 

(M  ^ = - 1)  shows  the  phase  velocity  which  this  family  of  solutions  begins  at 
(a  = 0)  as  a function  of  shock  Mach  number.  For  this  family  of  solutions 
to  have  amplified  solutions,  cr  < 1 and  for  neutral  solutions,  Cf  = 1 

From  Figure  6 can  be  seen  that  a supersonic  relative  flow 

region  for  subsonic  disturbances  with  c^  s 1 can  only  exist  when  M^  s 1.8 

in  shock  tube  boundary  layers.  Thus,  regular  amplified  subsonic  solutions 

are  not  possible  below  this  Mach  number.  Some  typical  regular  family 

solutions  are  shown  in  Figures  12  and  13  for  M = 1.  5,  2,  and  4.  In 

s 

Figure  12,  the  phase  velocity  (c^)  versus  wave  number  fa  ) is  shown  for  both 

viscous  and  inviscid  solutions.  As  shown  in  Figure  13,  the  viscous  solutions 

are  all  damped  whereas  the  inviscid  solutions  are  all  neutral.  The  only 

exception  is  at  M =4  where  the  inviscid  solutions  show  a very  small  amount 

of  amplification  for  a > 3.5.  For  a > 3.  5,  the  phase  velocity  (c^)  is  less 

than  one.  According  to  Mack,  a necessary  condition  for  an  amplified  inviscid 

regular  solution  to  exist  is  that  c^_  be  less  than  1.  Calculations  performed 

to  date  indicate  that  an  amplified  regular  viscous  solution  can  likewise  only 

exist  if  c <1.  In  Mack's  terminology,  the  solution  for  c = 1 is  called  the 
r r 

regular  neutral  solution  and  the  corresponding  value  of  a is  given  the  symbol 

a ir 
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STREAMWISE  WAVE  NO  O, 

Figure  10.  Effect  of  Wave  Angle  (i fi ) on  Amplification  (cj. 
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Figure  11.  Effect  of  Wave  Angle  (0)  on  Phase  Velocity  (c^). 
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The  viscous  counterpart  to  a j , denoted  by  30  lar§e 

that  extremely  small  computational  step  sizes  would  be  needed  to  treat  a 

condition  of  c si  for  the  regular  viscous  solution.  Additionally,  the 
r 

frequency  associated  with  this  solution  would  be  in  the  megahertz  range 
which  is  outside  the  range  of  frequencies  that  are  usually  considered  to  be 
of  importance  in  transition. 

With  one  exception  (to  be  discussed  in  the  following  paragraphs), 
all  amplified  solutions  available  in  the  literature  for  laminar  boundary  layer 
flows  are  encompassed  within  the  types  of  solutions  discussed  in  the  previous 
paragraphs. 

3.  4.  2 Supersonic  Outgoing  Solutions 

Mack  has  shown  in  Reference  6 that  with  sufficient  wall  cooling, 
amplified,  supersonic,  outgoing  solutions  belonging  to  the  c^  family  of 
solutions  can  exist.  However,  these  solutions  have  wave  numbers  which 
are  greater  than  so  that  again,  high  frequency  solutions  are  the  only 
possible  amplified  solutions.  Thus,  no  attempt  was  made  to  determine  the 
stability  characteristics  of  these  types  of  solutions  since  they  appear  to  be 
important  only  for  hypersonic  Mach  numbers  (see  Figure  11.25  of  Ref.  6). 

3.4.3  Summary  of  Results  for  Outgoing  Solutions 

The  results  of  stability  computations  for  outgoing  disturbances 

in  shock  tube  boundary  layers  encompassing;  two  dimensional  and  three 

dimensional  disturbances,  subsonic  and  supersonic  disturbances,  and  multiple 

mode  solutions  of  both  subsonic  and  regular  families,  have  not  yielded 

solutions  which  can  explain  or  correlate  with  boundary  layer  transition 

measurements.  The  outgoing  solutions  represent  disturbances  which  originate 

within  the  boundary  layer.  Therefore,  it  is  concluded  that,  within  the  context 

of  a small  disturbance  linear  boundary  stability  theory,  outgoing  disturbances 

in  shock  tube  boundary  layers  are  completely  damped  out  for  the  range  of 

Mach  numbers  considered  in  this  study  (M  = 1.  5 to  5.0).  It  is  further  con- 

s 

eluded,  subject  to  the  limitations  on  the  validity  of  the  linear  theory,  that 
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extremely  high  transition  Reynolds  numbers  would  be  attainable  if  it  were 
possible  to  eliminate  all  environmental  disturbances  in  a shock  tube.  Even 
at  a shock  Mach  number  as  low  as  1.  5,  linear  stability  theory  predicts  a 

g 

critical  Reynolds  number  as  high  as  approximately  1.7x10  . 

3.  4.  4 Incoming  Disturbances 

A new  type  of  solution  to  the  stability  equations  was  discovered 
during  the  course  of  the  outgoing  wave  study.  These  solutions  show  some 
promise  of  providing  a mechanism  for  explaining  transition.  These 
new  solutions  are  incoming  wave  eigensolutions  and  were  originally  dis- 
covered when  it  was  noted  that  certain  damped  supersonic  incoming  wave 

eigensolutions  end  abruptly  in  a c versus  a plot  unless  one  allows  the  sign 

r 

of  the  real  part  of  the  characteristic  value  of  the  viscous  solution(which  is  a 
counterpart  to  the  inviscid  free  stream)to  switch  from  negative  to  positive. 

Before  discussing  amplified  incoming  supersonic  disturbances, 
a number  of  related  points  should  be  made.  First,  Mack  had  never  found 
any  incoming  neutral  supersonic  solutions.  Second,  he  indicated  that  if 
they  did  exist,  they  would  not  have  the  importance  of  the  outgoing  neutral 
supersonic  solutions  (page  11-45  of  Reference  6).  Third,  the  amplitude 
of  an  incoming  amplified  supersonic  disturbance  would  have  to  increase 
with  increasing  distance  from  the  wall  (page  10-34  of  Reference  6).  The 
third  point  posed  some  conceptual  difficulties  at  first  since  only  those 
solutions  to  the  linear  stability  equations  which  are  at  least  bounded  at 
infinity  are  considered  to  be  valid.  This  conceptual  difficulty  was  resolved 
by  recognizing  that  the  boundary  layer  in  a shock  tube  does  not  extend  to 
infinity.  In  particular,  surface  roughness  or  turbulence  on  one  sidewall 
in  a shock  tube,  can  influence  the  boundary  layer  on  the  opposite  sidewall. 
Thus  solutions  which  have  amplitudes  increasing  with  distance  from  the 
wall  are  permissible  so  long  as  the  rate  of  growth  it  not  excessive.  This 

condition  is  generally  satisfied  by  the  viscous  counterpart  to  the  inviscid 

I 

solution  of  the  stability  equations.  The  initial  search  for  incoming  amplified 
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disturbances  was  prompted  by  an  unexpected  finding  made  during  investigation 

of  the  two  dimension  subsonic  disturbance  calculations  for  low  Reynolds  numbers 

and  low  wave  numbers  at  M = 3.0.  It  was  observed  that  as  cr  was  decreased 

s 

from  about  0.  3 to  0.  1 for  a 2-D  wave,  the  damping  rate  started  to  sharply 
decrease.  However,  3-D  waves  > 30°)  show  a gradual  decrease  in  damping 
and  then  a sharp  increase  as  a approaches  zero  as  illustrated  in  Figure  10. 

The  phase  velocity  of  the  2-D  subsonic  disturbances  shows  a steady  decrease 
first  approaching  and  then  dropping  below  c^-  1 - l/M^.  In  contrast  the  3-D, 

30°  wave  showed  a strong  increase  in  as  a approaches  zero  (as  shown  in 
Figure  11).  It  was  observed  that  the  character  of  the  2-D  disturbances  changed 
from  out-going  to  incoming  as  c^  dropped  below  1 - l/M^  . Decreasing  the 
Reynolds  number  from  R = 2500  to  1500  showed  that,  as  a dropped  below  0.075, 
the  damping  at  R - 1 500  became  less  than  the  damping  for  R = 2500.  These 
results  provided  the  first  clues  that  amplified  incoming  disturbances  might 
exist.  As  a was  decreased  to  about  0.05  and  the  Reynolds  number  decreased 
to  1000  it  was  observed  that  the  real  part  of  the  eigenvalue  corresponding 
to  the  inviscid  solution  had  to  be  allowed  to  change  sign  from  negative  to 
positive  as  it  passed  through  zero  in  order  to  maintain  an  incoming  distur- 
bance. A further  decrease  in  cy  at  R = 1000  yielded  the  first  amplified 
incoming  supersonic  disturbances  as  shown  in  Figure  14  and  15.  The 
results  presented  in  Figure  14  clearly  show  that  viscosity  is  destabilizing; 
that  is  the  amplification  factor  decreases  as  the  Reynolds  number  is 
increased.  Thus  these  incoming  solutions  cannot  be  found  from  an  inviscid 

stability  analysis.  Similar  calculations  were  performed  for  M =3.25 

s 

and  3.  5.  These  results  are  shown  in  Figures  16,  17,  18,  and  19. 

As  n is  decreased  at  constant  R,  in  Figure  16,  the  real  part 
of  the  characteristic  value  of  the  viscous  counterpart  to  the  inviscid  solution 
changes  sign  from  negative  to  positive  at  the  points  labeled  as  a "transitional" 
region.  For  values  of  & greater  than  a at  this  "transitional  region",  the 
amplitude  of  the  waves  decreases  as  y increses;  whereas,  for  rj  values  less 
than  ry  at  this  "transitional  region",  the  amplitude  of  the  waves  increases 
as  y increases.  These  new  incoming  wave  eigensolutions  must  be  interpreted 


Figure  14.  Effect  of  Streamwise  Wave  No.  (a)  on  Amplification  Factor 

(c. ) for  Varying  Reynolds  No.  (R  =Y R ) (106  Integration  Steps) 
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Figure  15.  Phase  Velocity  (c  ) Vs.  Streamwise  Wave  No.  (&  ) at  Varying 
Reynolds  No.  (R  =VR^  ) . 106  Integration  Steps,  = 3.0 


differently  from  outgoing  wave  eigensolutions.  For  che  incoming  wave 
eigensolutions,  positive  c.  means  that  the  strength  of  the  free-stream 
disturbance  is  growing  in  time,  while  negative  c^  implies  that  the  strength  of 
the  free-stream  disturbance  is  decreasing  in  time.  Thus,  c^  is  an  attribute 
of  the  free-stream  disturbance  and  is  not  related  by  itself  to  the  instability 
of  the  boundary  layer. 

The  physical  significance  of  amplified,  supersonic,  incoming  wave 
eigensolutions  is  associated  with  the  special  characteristic  which  these  waves 
possess,  namely,  these  waves  do  not  produce  a reflected  wave  in  the  boundary 
layer.  That  is,  while  an  ordinary  (noneigensolution)  incoming  wave  requires 
the  presence  of  a reflected  wave  at  the  wall,  the  incoming  wave  eigensolutions 
do  not. 

The  important  physical  significance  of  these  incoming  wave  eigen- 
solutions is  related  to  the  energy  transfer  associated  with  the  incoming  wave. 
All  of  the  energy  transported  into  the  boundary  layer  through  the  incoming 
wave  remains  in  the  boundary  layer;  there  is  no  reflected  wave  to  carry  energy 
out  of  the  boundary  layer. 

There  are  finite  Reynolds  stresses  associated  with  these  incoming 
wave  eigensolutions.  Thus,  these  incoming  wave  eigensolutions  provide  a 
mechanism  for  energy  redistribution  in  the  boundary  layer.  Since  there  is 
a net  energy  transport  from  the  disturbance  to  the  mean  flow,  these  incoming 
wave  eigensolutions  offer  a potentially  important  mechanism  for  changing  the 
character  of  the  mean  flow  profile. 

These  incoming  wave  eigensolutions  have  been  found  to  exist  for 
both  types  of  supersonic  disturbances,  i.e.,  for  c^<  l-l/M^and  c^l+l/M^. 
For  incoming  wave  eigensolutions  with  c^<  l-l/M^,  an  especially  interesting 
phenomena  has  been  found.  Two  separate  families  of  solut'  .is  have  been 
found  to  exist.  One  family  has  its  origin  at  c^=  1-1 /M^  and  <*  = 0 while  the 
other  family  originates  from  subsonic  solutions.  Above  a certain  Reynolds 
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number  labeled  "changeover  Reynolds  number"  (R^),  the  two  families 

switch  their  relative  positions.  This  phenomena  is  shown  in  Figures  20  and 

21  for  M =3.  25.  In  Figure  22,  the  variation  of  R , is  shown  as  a function 
s ch 

of  shock  Mach  number.  The  physical  importance  of  this  changeover  Reynolds 

number  appears  to  be  that  for  Reynolds  numbers  less  than  the  changeover 

Reynolds  number  the  forced  response  analysis  shows  that  a relatively  strong 

reflected  wave  is  present,  i.  e.  , significant  energy  is  carried  out  of  the 

boundary  layer.  For  Reynolds  numbers  greater  than  R , the  strength  of  the 

reflected  wave  is  relatively  weak.  Thus,  when  the  Reynolds  number  is  greater 

than  R , the  boundary  layer  is  more  susceptible  to  mean  flow  profile  mod- 
ch 

ification  than  when  R < R , . 

eh- 

The  variation  of  R , versus  M shown  in  Figure  22  does  follow 
ch  s 

the  general  trend  of  the  first  transition  reversal  observed  by  Boison.  Thus, 
it  would  appear  that  this  "changeover  Reynolds  number"  might  have  the  same 
physical  significance  for  incoming  waves  as  the  critical  Reynolds  number 
has  for  outgoing  waves.  That  is,  the  "changeover  Reynolds  number"  seems 
to  represent  a Reynolds  number  which  must  be  reached  before  significant 
transition  mechanisms  can  take  effect. 

The  transition  measurements  shown  in  Figure  1 suggest  that  the 

changeover  Reynolds  number  associated  with  supersonic  and  sonic  waves 

having  c < l-l/M^  maybe  responsible  for  the  transition  reversal  which 

occurs  in  the  vicinity  of  T / T ~ . 35.  Some  other  mechanism  must  be 

we 

responsible  for  transition  for  T /T  >.  35.  One  possibility  is  that  there 

we 

may  be  a similar  changeover  Reynolds  number  phenomena  associated  with 

waves  having  c > 1 + 1/M-.  This  possibility  was  not  investigated, 
r “ L 

3.  5 FORCED  RESPONSE  ANALYSIS 

Neutral  incoming  wave  solutions  (noneigensolutions)  have  been  used 
in  conjunction  with  neutral  outgoing  wave  solutions  by  Mack  ^ ^ to  determine 
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Figure  22.  Effect  of  Shock  Mach  Number  on  "Changeover" 

. Reynolds  Number. 

the  forced  response  of  the  laminar  boundary  layer  to  free  stream  super- 
sonic waves.  A similar  approach  was  applied  to  shock  tube  boundary  layers  in 
the  present  work.  This  effort  was  very  limited  in  scope  and  cannot  be 
considered  complete  at  this  point.  To  date  no  evidence  has  been  found  that 
the  forced  response  analysis  by  itself  can  explain  transition  t.i  shock  tube 
boundary  layer  flows.  Wall  cooling  has  the  same  effect  for  shock  tube 
boundary  layer  flows  as  it  has  for  supersonic  wind  tunnel  flows;  that  is, 

\yall  cooling  tends  to  weaken  the  response  of  the  boundary  layer  to  forced 
oscillations.  The  forced  response  analysis  however, yields  a particularly 
interesting  result  which  may  have  an  important  bearing  on  transition  in 
shock  tube  flows. 

Mack  has  shown  that  the  peak  response  of  the  boundary  layer  to 
supersonic  incoming  waves  of  the  noneigensolution  type  occurs  very  near 
the  leading  edge  of  the  boundary  layer.  The  limited  results  obtained  in 
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the  present  effort  also  show  this  type  of  behavior.  The  response  of  the 
boundary  layer  is  expressed  in  terms  of  the  strength  of  the  reflected  wave 
produced  by  the  incoming  wave.  Thus,  the  forced  response  of  the  boundary 
layer  produces  a reflected  wave  which  grows  in  strength  over  the  initial 
part  of  the  boundary  layer.  In  a shock  tube,  this  reflected  wave  will 
impinge  on  the  opposite  wall  and  can  possibly  excite  an  incoming  wave 
eigensolution  at  the  opposite  wall.  If  this  situation  does  occur,  then 
significant  energy  can  be  deposited  in  the  opposite  wall  boundary  layer.  The 
intriguing  possibility  thus  exists  for  explaining  how  boundary  layer  transition 
can  be  influenced  by  the  diameter  of  the  shock  tube.  The  larger  the  shock  tube, 
the  longer  the  time  required  for  the  growing  reflected  wave  to  impinge 
on  the  opposite  side  wall  boundary  layer. 


The  preceeding  comments  are  purely  speculative  and  obviously  further 
investigation  is  needed  to  support  this  conjecture.  A very  recent  set  of 

(15) 

experiments  conducted  by  Golobic  ' definitely  shows  that  the  diameter  of 
the  shock  tube  does  play  an  important  role  on  boundary  layer  transition. 
Figure  23  shows  the  boundary  layer  transition  measurements  reported  in 
Reference  15.  The  work  reported  in  Reference  15  was  performed  in  order 
to  obtain  transition  data  for  the  wall  cooling  range  of  0.  255  < T^/T^  < 0.  39. 
The  purpose  of  these  experiments  was  to  verify  the  existence  of  Boison's 
transition  reversal  loop,  (see  Section  1.2).  As  can  be  seen  in  Figure  23, 
the  results  of  Golobic  show  a marked  deviation  from  Boison's;  that  is 
from  the  expected  results.  The  primary  difference  between  Golobic  s 
shock  tube  and  Boison's  was  the  diameter  of  the  tube.  Golobic's  was 
17  inches  in  diameter  and  Boison's  was  4 inches  in  diameter.  In  general, 
Boison's  experiments  were  conducted  with  a shock  tube  which  should  have 
provided  a less  hostile  environment  to  the  boundary  layer  than  Globic's 
since  more  precuations  were  taken  to  minimize  wall  roughness,  alignment 
of  the  shock  tube  sections,  etc.  But  since  Golobic's  transition  measurements 
were  higher  than  Boison's  over  the  same  wall  cooling  range  it  must  be 
concluded  that  the  difference  in  diameter  of  the  shock  tubes  had  an  important 
effect  on  the  differences  in  transition  in  these  two  shock  tubes. 
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Figure  23.  Transition  Reynolds  Number  Correlation  Showing 
Results  from  Golobic's  Measurements 

3.  6 EFFECT  OF  VARYING  WALL  COOLING  ON  NEUTRAL 
STABILITY  OF  SUPERSONIC  INCOMING  WAVES 

A brief  investigation  was  undertaken  to  determine  the  sensitivity  of 

the  stability  of  supersonic  incoming  waves  to  a change  in  wall  temperature. 

These  calculations  were  performed  for  two  shock  Mach  numbers,  M =2.5 

s 

and  Mg  = 2.75.  The  effect  of  changing  the  wall  temperature  on  neutral 

stability  (c.  = 0.  0)  is  shown  in  Figure  24  for  M = 2.  5.  The  effects  of  small 
1 s 

changes  in  wall  temperatures, such  as  occur  for  instance  as  the  wall  becomes 

heated  by  the  flow  behind  the  shock,  is  seen  to  be  very  small.  The  results 

for  Mg  = 2.75  (not  shown)  similarily  show  a relatively  small  sensitivity 

to  changes  in  T / T when  the  Mach  number  is  held  constant, 
w e 
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Skin  friction  drag,  caused  by  turbulence  in  compressible  boundary 
layers, comprises  a significant  portion  of  the  total  aerodynamic  drag  ex- 
perienced by  subsonic  aircraft.  This  drag  component  can  be  reduced  con- 
siderably if  the  boundary  layers  can  be  maintained  laminar.  Hence,  phen- 
omena affecting  boundary  layer  transition  in  subsonic  flows  are  of  major 
importance,  especially  if  they  point  the  way  toward  developing  technology 
for  controlling  transition  and  delaying  its  onset.  The  stability  analysis  for 
shock -induced  boundary  layers  described  in  Section  III  have  been  extended 
to  treat  boundary  layers  in  steady,  subsonic  flows  such  as  are  produced  in 
conventional  wind  tunnels  and  in  subsonic  flight.  A description  of  this  analysis 
and  its  results  follows. 


4.1  FORMULATION  OF  THE  PROBLEM 


Solution  of  the  linear  stability  equations  for  steady  flow  subsonic 
wind  tunnel  boundary  layers  is  similar  to  that  discussed  in  Section  II  for 
shock  tube  boundary  layers.  The  major  difference  is  in  the  shape  of  the 
mean  flow  profiles.  The  shock  tube  profiles  are  of  the  Rayleigh  type  and 
the  wind  tunnel  profiles  are  of  the  Blasius  type  (see  Figures  4 and  5). 

Computer  generation  of  these  cooled  steady  flow  mean  flow  profiles 
presents  no  significant  difficulties.  All  that  is  required  is  specification 


of  the  dimensionless  enthalpy  at  the  wall  (0  ) given  by; 


* * , * * 

0 = (h  - h )/  (h  - h ) 

w w e o e 
e 


where, 


h = enthalpy  at  the  wall, 
w 


1 

i 

i i 

. " 
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h = enthalpy  in  the  free- stream, 
e 

h*  = stagnation  enthalpy  of  the  free- stream, 
e 

Once  these  mean  flow  profiles  have  been  determined,  they  are  used  as  input 
to  the  linear  stability  computer  program. 

4.2  RESULTS 

Neutral  stability  solutions  were  obtained  for  two  dimensional  and 
three  dimensional  disturbances  in  M = 0.603,  and  M = 0.80,  and  M = 0.90 
steady  flow  laminar  boundary  layers.  Calculations  were  performed  for  the 
following  thermal  conditions: 


1. 

insulated 

surface , 

2. 

T /T  = 
w e 

0.824, 

3. 

T /T  = 
w e 

0.759, 

4. 

T /T  = 
w e 

0.620. 

In  the  present  study,  second  mode  eigensolutions  to  the  viscous 
linear  stability  equations  were  not  considered  since  the  cimensional 
frequencies  associated  with  these  solutions  are  in  the  megahertz  range. 
Such  high  frequencies  are  not  physically  significant  for  subsonic  flows. 

Only  neutral  stability  calculations  were  performed  since  these  provide  a 
good  indication  of  the  general  stability  of  the  viscous  boundary  layer. 

Both  2-D  and  3-D  disturbances  were  considered. 

4.2.1  Two  Dimensional  Disturbances 

Results  indicate  that  for  neutrally  stable  two-dimensional  subsonic 
disturbances,  surface  cooling  increases  the  "critical  Reynolds  number" 
(Rcrit^*  Typical  results  are  shown  in  Figures  25  and  26.  Figures  25  (a) 
and  26  (a)  show  the  variation  of  the  neutral  stability  wave  number  (<y)  with 


BLASI US  PROFILE 
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Figure  25.  Effect  of  Wall  Cooling  on  Neutral  Stability  Conditions  for  2-D 
Disturbances  at  M = 0.603;  a)  Wave  Number,  b)  Phase 
Velocity. 


Reynolds  number  for  Mach  numbers  of  M = 0.603  and  0.90  respectively. 

Figures  25  (b)  and  26  (b)  denote  the  variation  of  phase  velocity,  c , with 

r 

Reynolds  number  at  the  same  Mach  numbers.  For  either  case,  as  the 

ratio  of  T /T  is  decreased,  the  critical  Reynolds  number  increases 
we 


monotonically  from  approximately  R = 300  (insulated  surface)  to  approxi- 


mately R - 20,000  (T  /T  = 0.620).  This  behavior  was  found  to  be 

w e 

characteristic  of  this  Mach  number  range. 

Dimensional  frequencies  typical  of  the  first  mode  neutral  disturbances 
found  in  the  Mach  number  range  from  0.6  to  0.9  are  presented  in  Figure  27 
for  M = O.Q.  In  general,  the  boundary  condition  that  results  in  the  highest 
degree  of  neutral  instability,  is  an  insulated  surface.  This  also  results  in 


the  highest  neutral  disturbance  frequency  (on  the  order  of  200  kHz). 


Figure  27.  Effect  of  Wall  Cooling  on  the  Neutral  Stability 
F requency  for  2 - D Distrubance s at  M^:  0.9 
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An  increase  in  surface  cooling  (or  a decrease  in  the  temperature  ratio 

(T  / T ) decreases  the  maximum  neutral  disturbance  frequency.  Physically, 

w e 

this  suggests  that  it  is  possible  to  limit  the  frequency  spectrum  of  poten- 
tially unstable  disturbances  by  increasing  surface  cooling.  Insulated 
surface  boundary  layers  are  relatively  unstable,  since  the  spectrum  of 
disturbance  frequencies  that  are  potentially  unstable  is  large.  Highly 
cooled  boundary  layers  are  more  stable,  since  the  spectrum  of  potentially 
unstable  disturbances  is  reduced.  Therefore,  if  the  frequency  spectrum  of 
disturbances  is  large  and  random  ("white  noise"),  a cooled  boundary  layer 
will  be  more  stable  since  only  a small  frequency  band  is  potentially  unstable. 


4.2.2  Three  Dimensional  Disturbances 

The  behavior  of  three  dimensional  subsonic  disturbances  with 
surface  cooling  was  also  investigated.  These  results  are  presented  in 
Figures  28  to  30  for  the  typical  case  of  M = 0.9.  Neutral  stability  calcu- 
lations were  performed  for  each  thermal  condition  previously  described 
for  the  two  dimensional  case.  The  streamwise  wave  angle  of  the  disturbance 
was  varied  from  0 to  87.5  • As  subsonic  neutral  disturbances  become 
increasingly  oblique  (increasing  wave  angle),  the  stability  of  the  boundary 
layer  increases  (see  Figures  28-30). 

It  is  apparent  from  Figures  28-30  that  stability  produced  by 
surface  cooling  is  a minimum  for  two  dimensional  disturbances  (wave  angle 
of  0°).  Thus,  2-D  disturbances  should  be  the  most  important  in  transition 
from  laminar  to  turbulent  flow. 

An  exception  to  the  above  findings  occurs  for  a highly  cooled 
boundary  layer  at  M = 0.9  and  T /T^  = 0.620.  Figure  31  displays  the 
results  for  this  case.  The  most  unstable  disturbance  is  no  longer  two 
dimensional,  but  three  dimensional,  and  occurs  at  a streamwise  wave  angle 
in  the  vicinity  of  30°  to  45°.  This  result  is  significant  in  that  it  differs 
from  Mack's  findings^*.  Mack  proposed  that  three  dimensional  disturbances 
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Figure  28.  Effect  of  Streamwise  Wave  Angle  on  Neutral  Stability  Conditions 
for  3-D  Disturbances  for  Insulated  Wall  at  M = 0.90; 
a)  Wave  Number,  b)  Phase  Velocity. 
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Figure  29.  Effect  of  Streamwise  Wave  Angle  on  Neutral  Stability  Condition 
for  3-D  Disturbances  at  = 0.  90  and  Tw/Te  = 0.  824; 
a)  Wave  Number,  b)  Phase  Velocity. 
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Figure  30.  Effect  of  Streamwise  Wave  Angle  on  Neutral  Stability  Conditions 

for  3-D  Disturbances  at  M - 0.90  and  T /T  = 0.759; 

oo  V/  c 

a)  Wave  Number,  b)  Phase  Velocity. 
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become  the  predominantly  unstable  disturbances  for  cooled  boundary  layers 

for  Mach  numbers  of  about  0.7  or  0.8.  This  appears,  however,  not  to  be 

the  case.  At  M = 0.80  (not  shown)  two  dimensional  disturbances  were  found 

to  be  the  most  unstable  for  all  values  of  T /T  analyzed.  At  M = 0.9, 

we 

three  dimensional  disturbances,  which  were  only  slightly  oblique  (i.e.  wave 

angle  of  about  30°),  were  found  to  be  the  most  unstable  (and  only  at 

T /T  = 0.620). 
w e 

Figure  32  is  a typical  neutral  stability  frequency  versus  Reynolds 

number  distribution  for  oblique  disturbances  at  M = . 90  and  for  an 

00 

insulated  surface.  Here  the  limiting  condition  on  the  maximum  neutral 
stability  frequency  is  the  two  dimensional  disturbance.  For  disturbances 
which  are  increasingly  oblique,  the  frequency  spectrum  available  to  a 
potentially  unstable  disturbance  becomes  narrower.  This  behavior  is 


REYNOLDS  NO.  ~ R (R=/R„) 


Figure  32.  Effect  of  Wave  Angle  on  the  Neutral  Stability  Frequency  for 
Three  Dimensional  Neutral  Disturbances. 
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similar  to  the  effect  of  wall  cooling  on  two  dimensional  disturbances. 

Three  dimensional  disturbances  of  constant  frequency  react  to  changes  in 
wave  angle  in  a manner  analogous  to  the  reaction  of  two  dimensional  dis- 
turbances to  changes  in  surface  cooling.  The  results  tend  to  confirm  the 
hypotheses  of  other  investigators  such  as  Dunn  and  Lin  ^ and  Mack  . 
Most  of  the  work  reported  by  Mack  conerning  the  effects  of  surface  cooling 
on  stability  was  conducted  with  an  inviscid  linear  stability  theory  and 
supersonic  flow.  He  concluded  that  increasing  the  surface  cooling  in  the 
boundary  layer  serves  to  stabilize  2-D  and  3-D  first  mode  disturbances. 
Higher  mode  disturbances  are  increasingly  unstable  with  increasing 
surface  cooling. 

In  summary,  the  linear  stability  analysis  of  cooled  subsonic  wind 
tunnel  boundary  layers  indicates  that: 

1.  Wall  cooling  has  a pronounced  stabilizing  effect  on  two  and 
three  dimensional  first  mode  disturbances. 

2.  Two  dimensional  disturbances  are  more  unstable  than 

three  dimensional  disturbances  for  Mach  numbers  of  M = 0.90 

or  less,  and  surface  cooling  ratios  greater  than  approximately 

T /T  = 0.62. 
w e 

3.  The  effect  of  Mach  number  on  two  and  three  dimensional 
disturbances  is  negligible  l°r  0.60  < M < 0.  90. 

4.3  COMPARISON  TO  SHOCK  INDUCED  FLOW  RESULTS 

For  purposes  of  general  comparison  between  the  stability  of  cooled 

subsonic  wind  tunnel  boundary  layers,  and  shock  induced  subsonic 

boundary  layers,  two  particular  cases  were  considered: 

1)  M = 0.603,  wind  tunnel  flow  and,  2)  M =1.5,  M = 0.603,  shock  tube 
oo  sc 

tube  flow.  The  behavior  of  both  cases  was  investigated  over  a range  of 
wall  cooling  conditions.  Results  of  these  calculations  for  points  of  neutral 


stability  are  presented  in  Figure  33.  The  subsonic  shock  induced  boundary 

layer  is  more  stable,  possessing  a higher  critical  Reynolds  number  for  a 

given  wall  cooling  ratio  than  its  wind  tunnel  counterpart.  This  difference 

is  attributable  to  the  differences  in  mean  flow  profiles.  The  shock  tube 

profile  has  higher  gradients,  and  hence  transfers  more  energy  out  of  the 

boundary  layer  than  the  steady-flow  Blasius  type  profile.  (This  phenomena 

(7) 

was  also  observed  by  Ostrach  and  Thornton'  \ ) Stability  is  influenced  not 
only  by  boundary  layer  cooling,  but  to  an  even  greater  extent,  by  the 
gradients  of  the  mean  flow  profile.  Thus,  direct  correlation  between 
experimental  transition  data  obtained  in  wind  tunnels  and  that  obtained  in 
shock  tubes  is  not  expected. 


Figure  33.  Comparison  of  the  Relative  Stability  of  Shock  Tube  and  Wind 
Tunnel  Boundary  Layers,  as  a Function  of  Critical  Reynolds 
Number  and  Wall  Cooling  Ratio. 
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4.4  EFFECT  OF  WALL  COOLING  ON  LOCAL  SKIN 
FRICTION  COEFFICIENT 

A comparison  between  laminar  and  turbulent  skin  friction  values  for 
the  specific  flow  parameters  under  consideration,  is  presented  in  Table  I. 
The  following  definitions  apply: 


w 


dU 

drj 


r?=o 


Dimensionless  viscosity  at  the  surface 
Dimensionless  velocity  gradient  at  the  surface 


C 

*L 

Cf  fRk 
T 


Laminar  skin  friction  parameter 
Approximate  turbulent  skin  friction  parameter 


It  is  apparent  that,  although  laminar  skin  friction  increases  with  wall  cool- 
ing, it  is  still  an  order  of  magnitude  less  than  the  turbulent  skin  friction  at 
the  same  Reynolds  number.  A potential  thus  exists  to  significantly  reduce 
aerodynamic  drag  on  aircraft  surfaces  by  providing  sufficient  cooling  to 
delay  the  onset  of  transition. 
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SECTION  V 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 
5.  1 SUMMARY 

The  work  described  in  this  report  ha^  showtf-that  wall  cooling  can  be 
an  effective  means  of  increasing  the  stability  of  laminar  boundary  layer  flows. 
The  stability  of  two  general  classes  of  cooled  boundary  layer  flows  was  in- 
vestigated; shock  tube  induced  boundary  layers,  and  cooled  subsonic  steady 
flow  boundary  layers  (such  as  that  on  a cooled  aerodynamic  surface  in  free 
flight  or  in  a wind  tunnel).  ^ 

A distinguishing  feature  of  wall  boundary  layers  in  shock  tubes  is 
that  the  wall  temperature  is  less  than  the  temperature  at  the  edge  of  the 
boundary  layer  for  all  shock  Mach  numbers.  In  this  report,  shock  Mach 
numbers  ranging  from  1.5  to  5.0  were  investigated.  The  results  showed 
that  the  cooling  of  the  boundary  layer  in  this  type  of  flow  is  sufficient  to 
supress  the  ordinary  types  of  disturbances  (2-D  and  3-D  subsonic  disturbances) 
considered  in  boundary  layer  stability  theory.  Higher  mode  solutions  to  the 
stability  equations  of  the  type  first  discovered  by  Mack  were  found  to  exist  in 
shock  tube  boundary  layer  flows.  The  frequencies  associated  with  these  higher 
mode  solutions  were  found  to  be  so  high  (in  the  megahertz  range)  that  they 
were  considered  to  be  unimportant  in  establishing  the  stability  of  shock  tube 
boundary  layers. 

The  primary  objective  of  the  reported  effort  was  to  determine  whether 

or  not  three-dimensional  and/or  higher  mode  solutions  correlate  with  trans- 

(4) 

ition  measurements  in  shock  tubes,  particularly  those  of  Boison  . When 
it  was  determined  that  these  types  of  solutions  would  not  correlate  with 
transition  measurements,  a search  f(  r other  types  of  solutions  which  might  show 
correlation  was  undertaken.  Two  new  classes  of  disturbance  solutions  were 
discovered  which  are  unstable  in  highly  cooled  boundary  layers.  Both  classes 
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originate  in  the  free  stream  in  the  form  of  supersonic  disturbances.  One 

class  corresponds  to  supersonic  disturbances  which  travel  slower  than  the 

free  stream  (c  < 1-1/M  ) and  the  other  class  corresponds  to  disturbances 
r co 

which  travel  faster  than  the  free  stream  (c  >H1/M  ).  Since  both  classes 

r oo 

of  disturbances  originate  in  the  free  stream,  these  unstable  disturbances 

are  not  important  for  determining  the  stability  of  cooled  boundary  layers  in 

a disturbance-free  environment  such  as  free  flight.  Conversely,  these  types 

of  disturbances  are  of  fundamental  importance  to  the  stability  of  shock  tube 

boundary  layers,  since  free  stream  disturbances  always  exist  in  shock  tube 

flows.  Free  stream  disturbances  also  exist  in  all  wind  tunnel  flows  so  that 

supersonic  incoming  disturbances  are  also  of  importance  in  determining  the 

stability  of  boundary  layers  on  models  in  wind  tunnels.  This  factor  is  especially 

important  for  cooled  models  when  the  cooling  is  sufficient  to  suppress  the 

subsonic  disturbances  which  would  otherwise  control  boundary  layer  stability. 

The  wall  cooling  required  to  suppress  subsonic  disturbances  was 

determined  for  steady  subsonic  flow  over  flat  plates  with  no  pressure  gradient. 

The  Mach  number  range  investigated  was  0.603  <M  <0.9  and  the  wall  tern- 

00 

peratures  range  extended  from  insulated  wall  conditions  to  T =0.62  T . The 

w e 

increase  in  critical  Reynolds  number  with  decreasing  T / T is  dramatic, 

w e 

as  is  shown  in  Figure  33.  Thus,  even  a moderate  amount  of  wall  cooling 

should  lead  to  a significant  increase  in  transition  Reynolds  number.  It  was 

also  found  that  two-dimensional  disturbances  are  more  unstable  (have  a lower 

critical  Reynolds  number)  than  three-dimensional  disturbances  over  this 

range  of  M and  T /T  , except  at  M =0.9  and  T /T  =0.62  where  an  oblique 
oo  we  oo  we 

wave  having  a wave  angle  of  about  30°  was  found  to  have  the  lowest  critical 
Reynolds  number.  In  addition,  it  was  shown  in  Section  4.4,  that  wall  cooling 
does  not  significantly  increase  the  local  skin  friction  coefficient  for  laminar 
fl  ow. 
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5.  2 CONCLUSIONS 

The  conclusions  cover  a wider  range  of  information  than  that  contained 
in  this  report  in  order  to  obtain  the  widest  possible  perspective  on  the  effect 
of  wall  cooling  on  boundary  layer  stability.  The  conclusions  do  not;  however, 
extend  beyond  the  range  of  flow  parameters  considered  in  this  report. 

The  first  five  conclusions  are  concerned  with  the  effect  of  wall 
cooling  on  subsonic  disturbances: 

1.  With  little  or  no  wall  cooling,  subsonic  disturbances  are 
the  principal  cause  of  transition. 

2.  Surface  cooling  stabilizes  subsonic  disturbances  and  delays 

transition  over  all  Mach  numbers  investigated  (M  < 2). 

00 

3.  Subsonic  disturbances  cause  transition  on  aerodynamic  surfaces 
in  wind  tunnel  flows  and  free  flight  unless  sufficient  surface 
cooling  is  imposed. 

4.  Shock  tube  flows  always  have  significant  wall  cooling. 

5.  Amplification  of  subsonic  disturbances  in  shock  tube  wall 
boundary  layer  flows  is  delayed  well  beyond  the  location  of 
measured  transition. 

The  wall  cooling  associated  with  shock  tube  flows  was  found  to  be  so 
effective  in  suppressing  both  two-dimensional  and  three-dimensional  dis- 
turbances that  these  types  of  disturbances  did  not  provide  any  correlation 
between  boundary  layer  instability  and  measured  transition  Reynolds  numbers. 
Therefore,  a search  for  other  types  of  unstable  waves,  which  could  account 
for  earlier  transition,  was  undertaken.  The  investigation  identified  a class 
of  supersonic  disturbances. 

From  a study  of  the  behavior  of  these  disturbances  it  was  concluded 

that; 

6.  Unstable  outgoing  supersonic  disturbances  exist,  but  their 
characteristic  frequencies  are  too  high  to  affect  boundary 
layer  stability. 


iJ 
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7.  Incoming  supersonic  disturbances  always  exist  and  have 
unstable  regions  over  all  flow  conditions  investigated. 

8.  Unstable  incoming  supersonic  disturbances  control  trans- 
ition in  shock  tube  flows. 

5.3  RECOMMENDATIONS 

The  results  in  this  report  show  that  wall  cooling  in  subsonic  and 
low  supersonic  Mach  number  (M<2)  has  a significant  potential  for  reducing 
skin  friction  drag  on  flight  vehicles  by  delaying  transition  from  laminar  to 
turbulent  flow.  The  results  of  this  investigation  were  obtained  for  flows 
with  no  pressure  gradient.  Obviously,  it  will  be  necessary  to  determine  if 
wall  cooling  can  prevent  laminar  flow  separation  and  transition  in  regions 
of  adverse  pressure  gradient  before  an  estimate  can  be  made  of  the  total 
aerodynamic  drag  reduction  which  can  be  achieved  through  wall  cooling 
on  flight  vehicles  (skin  friction  drag  + pressure  drag). 

The  following  program  for  further  research  will  be  necessary  to 
determine  the  overall  feasibility  of  boundary  layer  control  through  wall 
cooling. 

1.  A Wind  Tunnel  Investigation  of  Boundary  Layer  Transition 
on  Cooled  Flat  Plates  with  no  pressure  Gradient  for  sub- 
sonic Mach  numbers. 

2.  A Linear  Stability  Theory  Investigation  of  Subsonic  Laminar 
Boundary  Layers  on  Cooled  Surfaces  with  Adverse  Pressure 
Gradient. 

3.  An  Experimental  Investigation  of  The  Effect  of  Wall  Cooling 
on  Lift  and  Drag  on  Airfoils  at  Subsonic  Conditions.  This 
program  should  include  both  wind  tunnel  and  flight  testing. 

4.  A Systems  Analysis  to  Determine  the  Overall  Payoff  Which  Can 
be  Achieved  by  Integrating  Cooled  Aerodynamic  Surfaces, 
Supporting  Hardware,  and  the  Use  of  Cryogenic  Fuels  Such  as 
Hydrogen  into  a Flight  Vehicle. 
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APPENDIX  A 
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i.  * i 


CRITIQUE  OF  OSTRACH- THORNTON  ANALYSIS 


The  mean  flow  field  in  a shock  tube  is  basically  an  unsteady  flow. 

In  an  experiment,  the  instrumentation  is  fixed  in  the  laboratory  coordinate 
system  and  so  one  wants  to  understand  the  flow  from  the  point  of  view  of 
an  observer  fixed  in  space.  However,  from  the  point  of  view  of  an  observer 
riding  on  the  shock  wave,  the  flow  is  steady  (neglecting  small  changes  in  the 
free  stream  conditions  of  the  shocked  gas  due  to  the  presence  of  a growing 
boundary  layer  and  attenuation  of  the  shock  speed  itself).  Thus,  the  analytical 
formulation  of  both  the  boundary  layer  equations  for  the  mean  flow  and  the 
small  disturbance  equations  is  simpler  when  done  in  the  moving  coordinate 
system  where  the  mean  flow  is  constant  in  time  and  all  of  the  unsteadiness 
of  the  flow  is  due  to  disturbances.  The  correlation  of  boundary  layer 
stability  and  transition  data  with  analytical  predictions  however  is  easier 
done  in  the  fixed  coordinate  system  since  all  available  data  are  obtained 
with  sensors  fixed  in  the  laboratory  coordinate  system.  In  the  present 
analysis  the  mean  flow  is  solved  in  the  shock-fixed  coordinate  system 
and  then  transformed  back  into  the  laboratory-fixed  coordinate  system. 

The  entire  stability  analysis  is  formulated  in  the  laboratory-fixed 
coordinate  system. 


(7) 


In  the  Ostrach-Thornton  investigation  the  entire  analyis  was  carried 
out  in  the  shocked-fixed  coordinate  system.  The  small  disturbance  equations 


were  first  nondimensionalized  with  free,  stream-quantities  of  the  shocked  gas 


# sjs 

(U  , T , etc.  ).  Then,  basically  in  the  interests  of  expediency,  they 
e e 


normalized  the  mean  flow  velocity  profile  and  phase  velocity  in  a manner  so 
as  to  reduce  the  problem  to  one  in  which  the  normalized  velocity  profile  had 
a value  of  zero  at  the  wall  and  1 in  the  free  stream  corresponding  to  the 
classical  stability  analyses.  They  defined  new  quantities,  Q,  o'*  and  c 
as  follows 
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where  capital  letters  refer  to  dimensional  quantities  and  lower  case  letters 
refer  to  quantities  nondimensionalized  by  and  ( ) 1 signifies  differentiation 

with  respect  to  q , the  Blasius  similarity  variable.  Note  that  the  operation 
of  dividing  velocities  by  (1  - aj  when  these  velocities  were  originally  non- 
dimensionalized by  U*  is  equivalent  to  nondimensionalizing  the  velocities 
by  U*  - IT  in  the  first  place. 

The  disturbance  velocities  in  the  Ostrach-Thornton  analysis  were 
originally  nondimensionalized  by  U*  and  then  subsequently  divided  by  1 - u^. 
As  an  example,  the  normalized  value  of  the  amplitude  of  the  velocity 
fluctuation  in  the  streamwise  direction  f {their  notation)  became 


* 

f 


f/U 


(1-u  ) 
' w 


(u*  - tr ) 

' e w 


(A -2) 


where  here  f represents  a dimensional  quantity. 

Referring  to  Equation  (A-l)  above,  it  should  be  noted  that  the 

operations  performed  by  Ostrach  and  Thornton  as  indicated  in  Equation  (A-l) 

are  really  equivalent  to  transforming  the  mean  flow  velocity  and  phase 

velocity  back  into  a lab-fixed  coordinate  system.  This  transformation  is 

nerformed  when  U*  is  subtracted  from  the  mean  flow  velocity.  Then  they 
r w 

simply  nondimensionalized  all  velocities  by  the  velocity  of  the  fluid 
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relative  to  the  wall  (U^  - U^).  These  operations  clearly  lead  to  changes 

in  the  sign  of  the  quantities  ft  - c and  q'  compared  to  u - c and  u' 

respectively.  If  originally  u - c was  positive,  then,  after  the  operations 

of  subtracting  u individually  from  u and  c and  dividing  by  1 - u , ft  - c" 

W ° J yy 

becomes  negative  and  vice  versa  if  u - c was  originally  negative.  The  sign 

of  fl'  is  also  changed  from  that  of  u'  when  u'  is  divided  by  1 - u . There- 

w 

fore,  while  the  operations  indicated  by  Equation  (A-l)  appear  to  be  simply 
a normalization  procedure  used  to  obtain  a normalized  velocity  profile  that 
corresponds  to  the  classical  one  ranging  from  0 at  the  wall  to  1 in  the  free- 
stream,  it  must  be  recognized  that  these  operations  imply  more  than  just 
a simple  normalization.  The  operations  indicated  in  Equations  (A-l)  and 
(A -2)  (corresponding  to  Equations  7a  and  7b  in  Reference  7)  really  corres- 
pond to  first  transforming  the  small  disturbance  equations  back  into  the 
laboratory  coordinate  system  and  then  nondimensionalizing  all  velocities 
( by  the  velocity  of  the  fluid  relative  to  the  wall.  Therefore,  in  evaluating 

the  Ostrach  and  Thornton  analysis,  it  must  be  kept  clearly  in  mind 
which  coordinate  system  each  of  the  quantities  related  ter  velocities  is 
defined  in.  Q and  c are  properly  defined  in  a laboratory  fixed  coordinate 
system  and  thus  c is  the  nondimensional  phase  velocity  relative  to  a 
stationary  observer.  That  c is  thus  defined  is  not  stated  or  immediately 
evident  in  Reference  7. 

J 

, ' It  is  very  easy  to  become  thoroughly  confused  when  reading 

Reference  7 since  all  quantities  related  to  velocities  seem  to  have  the  wrong 
sign.  Ostrach  and  Thornton  apparently  recognized  that  these  quantities 


did  not  have  the  correct  sign  for  they  invented  new  sign  conventions  in 

order  to  "make  things  coxne  out  right".  They  defined  a new  wave  number 

a = - a and  the  Reynolds  number  Re,,  was  defined  as  I 1 - u I Re  instead 

1 ' w 1 

of  (1  - u^)  Re.  They  do  not  state  it  but  in  their  analysis  the  wave  number 
n the  transverse  direction  (3  would  also  have  to  be  defined  with  a negative 
;n.  Had  they  considered  3-D  waves  instead  of  only  2-D  waves,  it  is  quite 
t-  le  that  noticeably  incorrect  results  would  have  resulted  and  steps 
• > rre«  t the  analysis. 

* L * 7 
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It  is  now  clear  to  us  that  the  Ostrach-Thornton  analysis  is  basically 
properly  posed  (within  the  accuracy  of  the  Dunn-Lin  approach)  if  one 
untangles  all  of  the  sign  inconsistencies.  One  then  finds  that  they  have 
really  stated  the  problem  (equations  la  through  5a  in  Reference  7)  for  the 
situation  pictured  in  the  following  sketch 


Y 


and  that  all  disturbances  are  left  running  waves.  In  the  above  sketch  then, 

.w  JU 

- 3 c it'  U . and  U would  all  be  negative  quantities  as  they,  of  course, 

-*•  ’ M ’ ’ e w 

should  be  for  this  physical  situation.  The  input  mean  flow  profile  would  then 

be  input  as  negative  numbers  for  u and  positive  numbers  for  u.  Then  upon 

nondimensionalizing  with  respect  to  U*  - U " (a  positive  number)  fi  and  c 

e W _ 

would  be  positive  and  q - c would  have  the  correct  sign.  c would  indeed 

~ It  apoears  that  O-S  used  a positive  velocity  profile  for  U!,s  With  U*’ 
negative  (corresponding  to  their  Figure  7 for  steady  flow  system)  but  con- 
sidered Uw  and  Ug  also  to  be  positive  quantities.  They  considered  c to  be 
positive,  but  again  U*\,  and  U*  were  considered  to  be  positive.  Thus,  their 
q - c and  f.’  had  the  incorrect  sign  for  our  sketch.  But.  then  they  redefine  their 
disturbance  velocity  amplitudes  to  have  negative  values  which  tended  to  undo 
the  incorrect  signs  on!)-c  and  Q1.  Finally,  after  they  invented  a negative 
sign  for  a (and  by  implication  negative  signs  for  Band  c),  their  equations  la 
through  5a  became  applicable  to  the  physical  situation  shown  above.  When  they 
wrote  their  final  equations  (Eqns.  9a  through  14  in  Reference  7),  they  had 
in  essence  retransfornied  their  disturbance  velocities  and  redefined  a and 
to  correspond  to  the  flow  situation  which  is  just  the  opposite  of  our  sketch, 
but  corresponds  to  their  Figure  7. 


*» 


be  the  phase  velocity  relative  to  an  observer  in  the  laboratory-fixed 
coordinate  system  corresponding  to  the  sketch  which  is  really  what  the 
Ostrach-Thornton  phase  velocities  c are. 


In  Figure  34,  the  movement  of  a given  particle  is  tracked  in  the 

two  different  coordinate  systems  corresponding  to  the  shock-fixed  coordinate 

system  for  which  the  Ostrach-Thornton  analysis  is  really  applicable.  This 

figure  is  particularly  useful  for  visualizing  the  basis  for  the  development 

of  the  Reynolds  number  based  on  velocity  of  a particle  relative  to  the  wall 

and  length  equal  to  the  total  distance  that  a particle  is  in  motion  relative  to 

the  wall  at  a given  time  t.  Note  that  for  a time  interval  At=t^-t^  after  the 

shock  has  oassed  a given  oosition  (x  ) the  particle  that  is  at  the  oosition 

g 

x has  been  in  motion  for  the  time  t,  and  has  covered  the  distance  equal  to 
g 2 

U2  t?.  Also,  note  that  (all  quantities  are  considered  to  be  dimensional  in 
what  follows  ) 


At  = 


Xs  2 " X*1 


Xs 2 ‘ X£ 


and  that  At  = 


XP2  * Xpl 


x - x 


where  t is  the  time  after  the  shock  moves  past  a given  point  on  the  shock 
tube  wall.  Equating  these  two  expressions  for  ^t  yields 


XS2  * X 


x * xa. 

g ?1 


VUs  ‘ U2> 


U2  A t 


(Us  - u2) 


At  = 


which  states  that  in  the  time  At  after  the  shock  passes  a given  point  on  the 

shock  tube  wall,  the  particle  at  position  x at  time  t has  been  in  motion  for 

S 2 

the  time  given  by 
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or  Re v = (U  - U )2  U At/U 
X s e s e 

which  is  the  same  as  the  Reynolds  number  based  on  the  total  distance  that  a 
particle  travels  and  its  velocity  in  the  laboratory-fixed  coordinate  system. 

Therefore,  it  appears  that  the  end  result  of  the  Ostrach-Thornton 
formulation  of  the  stability  equations  is  correct  within  the  accuracy  of 
the  Dunn-Lin  approach  if  one  interprets  their  phase  velocity  c as  the 
phase  velocity  in  a shock-fixed  coordinate  system.  It  is  very  clear  that 
their  formulation  of  the  stability  equation  does  lead  to  the  correct  form 
of  the  parameters,  Mach  number  and  Reynolds  number,  for  a stability 
analysis  of  the  boundary  layer  formed  behing  a moving  shock  wave. 
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APPENDIX  B 

DESCRIPTION  OF  THE  BOUNDARY  LAYER  STABILITY 
COMPUTER  PROGRAM 

B.  I STABILITY  EQUATIONS  FOR  THE  VISCOUS  STABILITY  COM- 
PUTER PROGRAM 

The  complete  parallel  flow  linear  stability  equations  used  to  generate 
the  results  presented  in  this  report  can  be  tound  in  Reference  16.  The  stability 
equations  presented  in  Reference  16  contain  terms  which  arise  because  of  non- 
parallelism of  a boundary  layer.  These  terms  were  not  included  in  the  version 
of  the  stability  computer  program  used  to  generate  the  results  presented  in 
this  report.  The  terms  included  in  Reference  16  which  are  not  included  for 
parallel  flow  investigation  are  as  follows: 

(a)  All  terms  involving  V (mean  flow  component  of  velocity 
in  the  y-direction)  and  all  derivatives  of  V. 

(b)  All  streamwise  derivatives  of  mean  flow  quantities. 


The  stability  equations  in  Reference  16  also  include  the  effect  of 
the  bulk  coefficient  of  viscosity,  denoted  by  5.  In  the  present  investigation 
{ was  taken  to  be  zero. 

The  complex  amplitude  functions  of  the  flow  variables  are  defined 


by 


(u',  v',w',  p* ) = (f.  Of  0 , h,  n ) exp  [i(o  j x + » 3 z -U't)] 


(B.  1) 


When  equations  of  the  type  expressed  by  Equation  B.  1 are  substituted  into  the 
linearized  small  disturbance  equations  for  parallel  flow, an  eighth  order  system 
of  linear  differential  equations  is  obtained.  The  eighth  order  system  is  setup 
as  follows: 


Eb..Z.  (i=  l,2,o. 

U 1 


8) 


(B.  2) 
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where 

z j = f,  z2  = f = z'1,  z3  = 0 , z4  = tt/m^,  z5  = e 

z6  = e'  = z*5,  z?  = h,  zg  = h*  = z'? 

and  where  the  row  index,  i,  in  Equation  (B.  2)  is  set  up  in  the  following 
order 


Row  Index 
1 

2 

3 

4 

5 

6 

7 

8 


Equation 
Z’l  = Z2 

First  momentum  (x-direction) 
Continuity 

Second  momentum  (y-direction) 


Energy 


Third  momentum  (z -direction) 


The  above  formulation  of  the  stability  equations  given  above  is  identical 

2 

to  Mack's  except  that  he  defines  Z^  = ir/y  M#  . His  system  of  stability 
equations  also  has  one  further  minor  difference  and  that  difference  involves 
the  definition  of  the  bulk  viscosity  coefficient.  Mack  uses  the  Lees -Lin 
definition  which  is  3/2  times  the  usually  defined  value, 

The  system  of  equations  (Equation  (1))  is  written  in  the  form 


where 


Z. 

l 


D..  Z. 
ij  J 


8 


D.. 

ij 


2 dikbkj 

Ir  - 1 


( i=l, . . . , 8) 


(B.  3) 


(B.  4) 
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For  the  parallel  flow  version  of  the  stability  equations,  the  D_  can 
be  obtained  analytically  whereas  for  the  non-parallel  flow  version  the  D.. 
must  be  obtained  numerically.  The  values  for  the  parallel  flow  stability 
equations  are  given  below. 


The  Z^  equation  has  one  non- zero  coefficient 


D(l.  2)  =1 


(B.  5) 


The  equation  has  six  non-zero  coefficients 

D(2,  l)  = (ui  ^ +a  3)  + 1 ai  R(U'C) 

M T 


(B.6) 


D(2»  2)  = - — 

|i  dT 


(B.7) 
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D(2-3)=  + + ‘ <b-8> 
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D(2  4)  = - ( ~ + k ) a“M  (U-c)  + 

* '3  1 » y n 


(B.  9) 


D(2,5) 


<U-c)(f )-U  "Xl£_U_iiL  {B,10) 


H dT2  ^ dT 


D(2'6)=-f  ^ 


(B.  11) 


The  Z^  equation  has  5 non-zero  coefficients 


D(3,l)  = - i 


(B. 12) 


D(3,  3)  = T'/T 


(B. 13) 


D(3,  4)  = - i (U-c)  M 


(B.  14) 


D(3,5)  = — (U-c) 

D(3,7)  = -i^3 


(B. 15) 
(B.  16) 
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The  Z4  equation  has  no  non- zero  coefficients.  The  coefficients  are 


written  in  terms  of  the  factor  L,  defined  as 


R 


■+  i ( T + * ) M2  (U-c) 

u y fj,  3 00 


(B.  17) 


The  coefficients  are 


D(4,l)  = 


- 1 


r + ( — + k ) ^-] 

1 dT  ' 3 ' T J 


(B. 18) 


D(4,2) 


(B. 19) 


1 2 2 / 4 (T'f  du  T" 

D(4,  3)  = ~[-(«i  +a3  )+(—+«  ) — 


R a. 


- 1 


(U-c)] 


(B.  20) 


D(4,  4)  = -r(4-+«>  [0(3,4)  (^  |^+^)-iU’Mf]  (B.  21) 


L ' 3 ’ ' ' L"'"’  ■'  1 ji  d T T 
= JL  r/  — + * \ — -^-D(3,5)+  iU  ( 


D(4»  5)-  L [(  3 +*  > ^ dT 


_1_  d_ 
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(B.  22) 


D(4.6)  --i-  (-f-+«  ) 


(B. 23) 


a . 


D(4,  7)  - D(4,  1 ) 


a 


(B.  24) 


a . 


D(4,  8)  = D (4,2) 


a . 


(B.  25) 
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equation  has  a single  non -zero  coefficient 


equation  has  5 non -zero  coefficients 
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D (8,4)  = - ( j-  + K > aia3(U-c)Mf  + i y y- 


(B.  34) 


d (8,5)  My+Odjttj  m^r) 


(B.  35) 


D (8,  7)  = ( a ^ + a 3 ) + i«, 


1 M 


(B. 36) 


D (8,8) 
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(B. 37) 


In  these  equations 


* = ?/  M 

and  is  a constant. 

The  boundary  conditions  are 


(B.  38) 


ZY  (0)  = 0,  Z3(0)=0,  Z5(0)=0,  Z?(0)=0 


(B. 39) 


and  as  tj  -»  » , the  amplitude  functions  must  be  bounded,  or  for  the  case 
incoming  supersonic  waves,  the  amplitude  functions  must  be  bounded  for 

D>>  n.« 


The  problem  as  formulated  above  is  such  that  the  boundary  conditions 
themselves  are  not  sufficient  to  establish  any  solution  of  the  eight  differential 
equations  other  than  the  trivial  solution  of  zero.  Therefore,  the  problem  is 
an  eigenvalue  problem,  that  is,  nonzero  solutions  which  satisfy  the  boundary 


conditions  exist  only  for  certain  combinations  of  the  parameters  or^  , , R,  c. 


The  general  method  used  to  determine  the  eigenvalues  and  eigensolutions  to 
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the  above  system  of  equations  is  given  in  the  next  section.  The  computer 
program  used  to  implement  the  solution  is  given  in  Section  B.  3. 

B.  2 METHOD  OF  SOLUTION 

The  difficulties  of  solving  the  boundary  layer  stability  equations  are 
well  known.  Chapter  III  of  Reference  17  can  provide  the  reader  with  an 
elementary  understanding  of  the  main  difficulty  involved  in  numerically 
integrating  the  system  of  equations  represented  by  Equation  B.  3.  The  main 
difficulty  is  associated  with  the  problem  of  controlling  parasitic  error  growth 
in  the  solution  of  two -point  boundary  value  problems.  In  the  BLSTAB  program, 
the  Gram-Schmidt  orthogonalization  procedure  is  used  to  overcome  this 
problem.  The  reader  is  referred  to  Reference  18  (pp.  98-103)  for  a concise 
demonstration  of  the  problem  of  parasitic  error  growth  and  how  it  can  be 
overcome  by  the  use  of  the  Gram-Schmidt  orthogonalization  procedure. 

Chapter  IV  of  Reference  19  also  demonstrates  the  use  of  the  Gram-Schmidt 
orthogonalization  procedure  for  solving  two -point  boundary  value  problems. 

The  details  of  the  use  of  the  Gram-Schmidt  orthogonalization  procedure  to 
determine  the  eigenvalues  and  eigensolutions  of  Equations  B.  3 have  been 
presented  in  Reference  11. 

The  method  used  to  determine  the  eigenvalues  of  Equations  B.  3 
is  outlined  as  follows.  Four  linearly  independent  solutions  to  Equations 
B.  3 are  generated.  These  four  independent  solutions  are  chosen  to  satisfy 

the  boundary  conditions  at  Tl>>  7)  . A linear  combination  of  the  four  in- 

0 

dependent  solutions  is  then  determined  such  that  the  velocity  amplitudes 
are  zero  at  T)  =0  and  the  pressure  amplitude  at  T|  =0  is  1.0.  Using  this  same 
linear  combination,  the  value  of  0(0)  is  evaluated.  If  0(0)  is  sufficiently 
close  to  zero,  the  values  of  o^,  and  c which  were  used  to  generate 

the  four  independent  solutions  constitute  ' .>  eigenvalue  set.  The  details  of 
the  methods  outlined  above  are  discussed  in  the  following  paragraphs. 


B.  2.  1 Free  Stream  Solution 


In  the  free  stream  U = T=fx=l  and  U,=U"  = T'  = T"  = 0.  All  of  the 
coefficients  in  Equations  B.  3 are  constant  and  only  the  following  coefficients 
are  non-zero:  D (2,  1),  D (2,4),  D (2,  5),  D (3,  1),  D (3,4),  D (3,  5),  D (3,  7), 

D (4,2),  D (4,3),  D (4,6),  D (4,8),  D (5,6),  D (6,4),  D (6,5),  D (8,4),  D (8,5), 
and  D (8,  7).  An  analytic  solution  to  Equation  B.  3 can  be  obtained  by  ele- 
mentary methods.  The  easiest  way  to  obtain  the  analytic  solution  is  to  re- 
write Equations  B.  3 as  a system  of  four  second  order  equations  in  terms  of 
Zj,  Z^,  Zj.,  and  Z^.  This  system  of  four  second  order  equations  has  the 
form 


w.  = £ b ..  W . (i=l...  ,4) 


(B. 40) 


iere  W =Z  , W =Z  , W =Z  , and  W =Z  and  the  non-zero  b::  are  given 

11  Z 4 3 5 4 7 -1 


in  terms  of  the  D by: 


bll=b44 


D (2,  1) 


(B.41) 


b12  = D (2,4) 


(B.  42) 


b j 3 = D (2,5)  (B.  43) 

b = D ( 2,  4)  x D (4,  2)+D  (3,4)xD  (4,3)  +D  (6,4)xD  (4,6) 

+D  (8,  4)  x D (4,8)  (B.  44) 


b = D (2,  5)  x D (4,  2)+D  (3,5)xD  (4,3)  +D  (6,5)xD  (4,6) 

+D  (8,5)XD  (4.8) 


(B. 45) 


b32  = D (6,4) 
b33  = D (6,5) 


(B. 46) 


(B.  47) 


b42  = D (8,4) 


(B. 48) 


b43  = D (8,  5) 


(B. 49) 


b44  = D (8,7) 


(B.  50) 


r 1 


- i 


The  characteristic  values  of  the  system  of  Equations  B.  40,  which 
are  the  same  as  the  characteristic  values  of  Equations  B.  3,  are  given  by 


X j=  - (bjj)  ^ = -[(oj  + )+ia  jR(l-c)]1^ 

(B.  51) 

X 2=X  1 

) 1/2 

(B. 52) 

X3.  - [>/2(b22+b33H[l/4(b22.b33)2  + b32b23]1/2 

1 

11/2 

} 

1 

(B. 53) 

*4.  - j>/2(b224b33Hl/4,b22-b33)2+b32b23]1/2  j 

(B.  54) 

x5= 

(B. 55) 

x6=  -x3 

(B.  56) 

II 

00 

II 

1 

>* 

(B.  57) 

In  order  to  study  incoming  supersonic  waves  it  is  necessary  to  determine 
which  of  the  above  characteristic  values  represent  the  viscous  counterparts 
to  the  inviscid  solution  as  given  by  Equations  13,  14,  and  15  on  page  24.  To 
this  end,  it  is  instructive  to  evaluate  b^  and  b^  for  the  special  case  of 
* = 0 and  rr  = 3/4.  From  Equation  B.45  it  can  be  shown  that 


t>23=  Oj  R(l-c)  (1-4/3  rr  - k a) 


(B. 58) 


and  from  Equation  B.  44  it  can  be  shown  that 

a M 2 (1-c)2  [ Y -(4/3+  * )a  (Y  - 1 )] 

’ =(»,  + <£)- 


22  ""1 


[1+  i Oj  yM  (1-c)  (4/3+  « )] 

IT 


For  the  case  of  * =0  and  nr  = 3/4,  Equation  B.  58  reduces  to 


b23=° 


and  from  Equation  B.  53 
>3--(b22,'/2 


(B. 59) 


(B. 60) 


(B.61) 


‘ • V 
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.*  * 


) 


and  from  Equation  B.  54 

*4--‘b33)l/2 


(B.  62) 


Combining  equations  B.61  and  B.  59  yields  the  following  expression  for 
X when  « =0  and  rr  = 3/4 

3 . v 


12  2 2 
_ _ o , M (1-c) 

,2  2 1 ® 

(o'  ,+»  7 

1 [1+io,  4/3  YM  (1-c)] 

1 00 


(B. 63) 


which  for  large  R reduces  to 


3 = -a  l i[l-M*  (1-c)2]  + 

R— *oo 


(B.  64) 


Thus,  X , and  X are  the  viscous  counterparts  to  the  inviscid 
3 6 

characteristic  solutions  (Note  that  Equations  13-15  on  page  24  are  valid  only 
for  2-D  waves  whereas  Equation  B.  64  is  a 3-D  wave  result).  From  Equations 
B.  62,  B.  47,  and  B.  30  it  is  easy  to  show  that  X 4 is  given  by 


7 ? 1/2 

X 4=  -[(a  ^ + a 3)  + 3/4  i R(l-c)] 


(B.  65) 


when  k =0  and  <j=  3/4. 


The  characteristic  vectors  corresponding  to  each  of  the  eight 

characteristic  values  can  most  easily  be  obtained  by  first  determining  the 

characteristic  vector  components  corresponding  to  Z^,  Z^,  Z^,  and  Z^ 

from  Equations  B.40-B.  50,  and  then  using  the  relations  Z2=Z'  , Z^  = ZX^, 

Z =Z*  and  the  continuity  equation  to  define  the  other  components.  Let 

B^  represent  the  characteristic  vector  corresponding  to  the  characteristic 

value  X ..  Then  from  Equations  B.40-B.  50  the  following  system  of  equations 
J (i) 

for  the  B J are  obtained: 


b12  B2(j’  + b13  B3<i>  =° 


(B. 66) 
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(b22  -M)B2l))+b23  B3<i’-° 

(B.  67) 

b32  B2(i’  + ">33  B3<i’  -° 

(B. 68) 

b42  B2li>  + ”43  B3U’  + (b44 

(B. 69) 

1/2 

For  the  first  characteristic  value,  X -(bjj)  = 

arbitrary  as  is  B Therefore,  B J'  ^ and  B ^ ^ 

m J/2  (1)  . 

-(b44)  * B!  18 
must  both  be  zero. 

Observe  that  \ is  a double  root  to  the  characteristic  determinant.  How- 
ever, since  only  one  (and  not  two)  of  the  Equations  B.  66-B.  69  is  a linear 
combination  of  the  other  equations  for  Xj=X  it  follows  that  even  though 
there  is  a double  root,  there  are  still  two  linearly  independent  characteristic 
vectors  corresponding  to  the  double  root.  Thus,  B and  B are  defined 
as  (1, 0,  0,  0)  while  B(2)  and  B(?'  are  defined  as  (0,  0,  0,  1).  The  remaining 
characteristic  vectors  are  given  by  the  following  relations: 


R W 

B3  - 

-D(6, 4) 

(B. 70) 

Bz  = 

D(6,5)-\2. 

(B. 71) 

R (j> 
B1  " 

D(2,  4)  x B,(j)-D(2,  5)  x D(6,4) 

z 

(B.  72) 

- D(2,  1) 

D(8,  4)  x B2(j)  + D(8,  5)  X B^ 

R (j) 

(B.  73) 

B4  ‘ 

X2.  - D(8,  7) 
J 

The  characteristic  vectors  for  the  eighth  order  system  A ^ are 

then  defined 

as  follows: 

A (j) 
A1 

II 

«— •. 

Ui. 

(B.  74) 

A «> 

= 

(B.  75) 

2 

A (j) 
A3 

= [ D(3,  1)  x Bj(j)  + D(3,  4)  x B2(j)  + D(3,5)  x B3<j) 

+ 

D(3,  7)  x B4(j)]/X2j 

(B.  76) 

i 


1 
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V * B2<i)/Xj 


(B. 77) 


A (j)  = B (j)/x  . 
5 3 ] 


(B. 78) 


A (j)  - B (j) 
A6  ' B3 


(B. 79) 


A?'j)  = B4(j)/\. 


(B. 80) 


A (j)  - B (j) 
A8  " B4 


(B. 8 1 ) 


The  analytic  free  stream  solutions  provide  the  initial  conditions 
for  four  linearly-independent  solutions.  The  characteristic  vectors  corres- 


ponding to  the  characteristic  values  having  negative  real  parts  (X  j,  X X y 


and  X 4)  are  then  integrated,  each  separately,  across  the  boundary  layer 


to  produce  four  linearly-independent  solutions.  Each  of  these  four  solutions 
approach  zero  as  T)-*«  . When  amplified  incoming  supersonic  waves  are 


being  considered,  \ ^ and  it's  characteristic  vector  must  be  used  instead  of 


X and  it's  characteristic  vector  in  order  that  the  pressure  amplitude  have 


the  characteristics  of  an  incoming  Mach  wave. 


B.  2.  2 Solution  Inside  the  Boundary  Layer 


With  initial  conditions  at  T)  = T]  specified  by  the  analytic  solutions, 

6 


the  Equations  3.  3 are  then  numerically  integrated  across  the  boundary  layer 


from  Tj  =T[  to  Ti  =0.  A standard  fourth  order  Runge-Kutta  technique  is  used 


to  carry  out  the  numerical  integration.  The  four  separate  integrations  are 
carried  out  in  parallel  across  a specified  number  of  integration  steps.  Then 
the  Gram-Schmidt  orthogonalization  procedure  is  performed;  both  on  the 
initial  conditions  as  well  as  on  the  current  solutions.  When  the  wall  (T)  =0) 
is  reached,  a final  Gram-Schmidt  orthogonalization  is  performed  and  then 
a linear  combination  of  the  solutions  is  obtained  satisfying  the  boundary  con- 


ditions at  the  wall  (Z ^(0)=Z^(0)- Z^=0  and  Z4(0)=  1.0).  0(0)  is  then  obtained 


where  Qis  given  by 


1 
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0=C,  Z.(1)  + C_  Z_(2)  +C,  Z.(3)  +C.  Z.(4)  (13.82) 

15  25  35  45 

where  the  C.  are  the  combining  coefficients  obtained  by  satisfying  the  wall 
boundary  conditions.  If  0(0)  is  zero  within  a specified  tolerance,  an  eigen- 
value has  been  found.  If  0(0)  ^ 0,  the  Newton-Raphson  search  produce  is 
used  to  find  a combination  of  o j,  Ri  and  c which  will  yield  an  eigenvalue. 

B.  3 LISTING  OF  THE  BLSTAB  COMPUTER  PROGRAM 

The  computer  program  which  implements  the  methods  discussed 
in  Section  B.  2 is  called  the  Boundary  Layer  STABility  Computer  P.ogram. 
The  input  data  for  the  program  are  as  follows.  The  data  on  the  first  card 
are: 

FIRST  CARD  (all  quantities  in  110  format) 

Columns  1-10.  NSTEPS--the  number  of  integration  steps. 

Columns  1 1 -20.TTSTEP--the  number  of  integration  steps  between 

Gram -Schmidt  orthogonalizations; 

Columns  21  -30.  ORDER --the  order  of  the  differential  equation 

system  which  is  eight  (8)  for  the  version 
of  the  program  described  in  this  report.  , 

Columns  3 1 -4(LNPNTDA--if  not  equal  to  zero,  the  input  mean 

flow  profiles  will  be  printed  out. 

MEAN  FLOW  PROFILE  INPUT  CARDS 

The  program  requires  that  the  mean  flow  quantities  U , U 1 , , T, 

T',  and  T"  be  specified  at  the  wall  and  at  equal  intervals  of  one-half  of  the 
integration  step  size  (H/2).  Thus,  the  mean  flow  profile  data  must  be 
available  at  2xNSTEPS+l  locations  equally  spaced  across  the  boundary 
layer  at  intervals  of  H/2.  Two  input  cards  are  used  to  prescribe  the  mean 
flow  profile  data  at  each  point.  The  first  card  contains  U,U',  and  U"  in 
a 3E14.  7 format  in  the  order  given.  The  second  card  in  each  set  contains 
T,  T',  and  T"  and  Y,  the  location  within  the  boundary  layer  at  which  the 
data  are  given.  T,T',  and  T"  are  in  a 3E14.  7 format.  Y is  in  columns 
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71-80  in  a F10.  3 format.  The  order  in  which  the  mean  flow  profile  data 
must  appear  is  in  ascending  order  of  Y,  that  is,  the  first  two  cards  give  the 
data  at  Y=0  (T)=0),  the  next  two  cards  give  the  data  at  H/2,  and  so  on,  with 


the  last  two  cards  giving  data  at  T|  =T>  - The  last  data  set  is  used  in  sub- 

5 


routine  EIGEN  to  determine  the  analytic  solution  for  the  free  stream.  There- 
fore, the  last  two  cards  must  show  U(T)  J=l.  U'(T]  ^)=U"(,n  ^ )=0  on  the  second 


last  card  and  T(T1  . )=1.  0,  T'CD  . ) = 0.  T"(T1  ) = 0,  and  Y=T1  on  the  last  card. 

0 B 6 o 


H,  the  integration  step  size  is  calculated  within  the  program  by  dividing 


T]  by  NSTEPS.  The  user  must  insure  that  the  number  of  integration  steps 

6 


divided  by  ITSTEP  is  an  integer. 

DATA  AND  PARAMETERS  REQUIRED  FOR  EACH  PROBLEM  SET 


A set  of  four  cards  is  required  for  each  problem. 
First  Card  of  a Problem  Set 


Columns  1-20.  EPS,  the  tolerance  of  0(0).  If  | 0(O)|  < EPS,  an 
eigenvalue  has  been  found.  E20.  1 format. 


Columns  21-30.  NOIT2,  the  number  of  iterations  of  the  Newton- 
Raphson  search  procedure  which  are  permitted 
to  find  an  eigenvalue;  110  format. 


Columns  3 1 -40.  NPRINT,  a diagnostic  output  parameter.  If  NPRINT 
^ 0,  the  results  of  intermediate  computations 
associated  with  obtaining  the  free  stream  solution 
are  printed  out  each  time  a trial  solution  is  carried 
out;  110  format. 


Columns  41-50.  NPRINT  1,  another  diagnostic  output  parameter 
to  be  used  in  conjunction  with  NPRINT.  If  only 
the  characteristic  values  of  the  free  stream  solution 
are  desired  for  each  individual  solution,  then 
NPRINT1  ^ 0 and  NPRINT=0.  110  format. 


Second  Card  of  a Problem  Set  (all  quantities  in  an  Fll.  5 format) 


Columns  1-11.  M,  Mach  number  of  the  free  stream.  For  shock 
induced  flows,  enter  the  Mach  number  of  the  flow 
behind  the  moving  shock  wave  in  a laboratory-fixed 
coordinate  system. 


. i 

-•  i 
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Columns  12-22.  GAMMA,  specific  heat  ratio  cp/cv  at  free  stream 
temperature  (assumed  to  be  constant  throughout 
the  flow). 


Columns  23-33.  SIGMA,  Prandtl  number  of  the  free  stream  (assumed 
to  be  constant  throughout  the  flow). 

Columns  34-44.  TINF,  temperature  of  the  free  stream  in  degrees 
Rankine. 


Columns  45-55.  MU2,  bulk  viscosity  (denoted  by  k in  Sections 
B.  1 and  B.  2). 


Third  Card  of  a Problem  Set  (all  quantities  in  an  Fll.  5 format) 

Columns  1-11.  (ALPHA1)  ; real  part  of  o 

**  1 

rf  nf  a< 

1 


r * '1 

Columns  12-22.  (ALPHA1).;  imaginary  part  of  o. 


Columns  23-33.  (ALPHAS)^;  real  part  of  o^ 

Columns  34-44.  (ALPHA3L;  imaginary  part  of 
Columns  45-55.  R,  Reynolds  number  (R^/r^) 

Columns  56-66.  (PHV)^;  real  part  of  phase  velocity,  c^ 

Columns  67-77.  (PHV).;  imaginary  part  of  phase  velocity,  c^ 


Fourth  Card  of  a Problem  Set  (all  quantities  in  an  110  format) 


Columns  1-10.  1X1  The  seven  numbers  on  the  third  card  are 

stored  in  an  array  called  X.  1X1  denotes 
one  of  the  two  parameters  which  will  be 
varied  in  carrying  out  the  Newton-Raphson 
search  procedure. 

Columns  11-20.  1X2  Denotes  the  second  of  the  two  parameters 

which  will  be  varied  in  carrying  out  the 
Newton-Raphson  search  procedure. 

Columns  21-30.  ITEST.  a ^ can  be  input  in  two  ways.  If  ITEST 

equals  1,  the  input  number  for  o ^is  con- 
sidered to  be  j.  If  ITEST  ^ 1,  then 

the  numbers  read  in  for  o,  are  considered 
to  be  o y The  option  of  ITEST=1  is  used 
to  study  waves  which  have  a constant 
oblique  angle. 

Columns  31-40.  INT.  If  INT^  0,  the  eigensolution  for  the  region 

inside  the  boundary  layer  will  be  computed 
and  printed  out  after  an  eigenvalue  has  been 
found. 
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Columns  41-50.  INTOBL.  If  INTOBL^  0,  the  eigensolution  for 

the  region  outside  of  the  boundary 
layer  will  be  computed  and  printed  out. 

Multiple  problem  sets  can  be  included  in  the  input  data.  This 
concludes  the  description  of  the  input  data  for  the  program.  The  listing 
of  the  program  is  given  in  the  next  section. 
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